RADC-TR-07-160 
Final  Technical  bport 
October  1907 


ON  ADAPTIVE  CELL-AVERAGING  CFAR 
RADAR  SIGNAL  DETECTION 


Syracuse  University 


Mourad  Barkaf  and  Fra  mod  K.  Varshnay 


mm  FOR  PUBLIC  RBIASZ  DISWBUTWU  UBURUTED 


DTIC 

ELECTEI 
JUN  0  71988 


ROME  AIR  DEVELOPMENT  CENTER 
Air  Fere*  Systems  Command 
Griffin  Air  Fores  Bara.  NY  13441-5700 


3  8 


00  6 


mmmm 


This  report  has  been  reviewed  by  the  RADC  Public  Affairs  Office  (PA)  and 
is  releasable  *-to  the  National  Technical  Information  Service  (NTIS) .  At  NTIS 
it  will  he  releasable  to  the  general  public.  Including  foreign  nations. 

RADC-TRi»*8  7-160  has  been  reviewed  and  is  approved  for  publication. 


APPROVED: 


VINCENT  C.  VANNICOLA 
Project  Engineer  " 


4- 


APPROVED: 


DAVID  J.  PRATT,  Col,  USAF 
Director  of  Surveillance 


FOR  THE  COMMANDER: 


JAMES  W.  HYDE,  III. 

Directorate  of  Plans  &  Programs 


If  your  address  has  changed  or  if  you  wish  to  be  removed  from  the  RADC  mailing 
list,  or  if  the  addressee  is  no  longer  employed  by  your  organization,  please 
notify  RADC  (OCTS)  Griffiss  AFB  NY  13441-5700.  This  will  assist  us  in  main¬ 
taining  a  current  mailing  list. 


Do  not  return  copies  of  this  report  unless  contractual  obligations  or  notice 
on  a  specific  document  requires  that  it  be  returned. 


UNCLASSIFIED 
fV  ■£'LA??iTfOTiSN-4 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704-0188 


la.  REPORT  SECURITY  CLASSIFICATION 

UNCLASSIFIED _ 

2a.  SECURITY  CLASSIFICATION  AUTHORITY 

N/A _ 

2b.  DECLASSIFICATION /DOWNGRADING  SCHEDULE 

Jtth _ 

4.  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 
N/A 


6a.  NAME  OF  PERFORMING  ORGANIZATION 
Syracuse  University 


6c.  ADDRESS  [City,  State,  and  ZIP  Code) 

Department  of  Electrical  and  Computer  Engr 
Link  Hall 

Syracuse  NY  13244-1240 

8a.  NAME  OF  FUNDING /SPONSORING  I  8b  OFFICE  SYME 

ORGANIZATION  (If  applicable ) 

Rome  Air  Development  Center  OCTS 

8c  ADDRESS  (City,  State,  and  ZIP  Code)  1 

Griffiss  AFB  NY  13441-5700 


lb.  RESTRICTIVE  MARKINGS 
N/A 

3  DISTRIBUTION/AVAILABILITY  OF  REPORT 

Approved  for  public  release;  distribution 
unlimited . 

S.  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 

RADC-TR-87-160 

7a.  NAME  OF  MONITORING  ORGANIZATION 
Rome  Air  Development  Center  (OCTS) 

7b  ADDRESS  (C/fy,  State,  and  ZIP  Code)  . . 

Griffiss  AFB  NY  13441-5700 


8b  OFFICE  SYMBOL  9  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 
(If  applicable) 

.  F30602-81-C-0169 


16b.  OFFICE  SYMBOL 
(If  applicable) 


10.  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM 
ELEMENT  NO 

61102F 


PROJECT 

NO 


I  WORK  UNIT 
ACCESSION  NO 


11.  TITLE  (Include  Security  Classification) 

ON  ADAPTIVE  CELL-AVERAGING  CFAR  RADAR  SIGNAL  DETECTION 


12.  PERSONAL  AUTHOR(S) 

Mourad  Barkat,  Pramod  K,  Varshney _ 

13a  TYPE  OF  REPORT  Il3b  TIME  COVERED  " 

Final  from  Tun  to  nPr  86 

16  SUPPLEMENTARY  NOTATION 


14.  OATE  OF  REPORT  (rear.  Month,  Day)  15  PAGE  COUNT 

_  October  1987  _  172 _ 


17. 

COSATI  CODES  ) 

FIELD 

GROUP 

SUB-GROUP 

09 

03 

_ U _ 

.  04  . 

18.  SUBJECT  TERMS  ( Continue  on  reverse  if  necessary  and  identify  by  block  number) 

Fusion  Detection 

Surveillance  Estimation 


In  radar  signal  detection,  the  problem  is  to  automatically  detect  a  target  in  a  nonstationary 
noise  and  clutter  background  while  maintaining  a  constant  probability  of  false-alarm. 
Classical  detection  using  a  matched  filter  receiver  and  a  fixed  threshold  is  not  applicable 
due  to  the  nonstationary  nature  of  the  background  noise.  Therefore,  adaptive  threshold 
techniques  are  needed  to  maintain  a  constant  false-alarm  rate  (CFAR) .  One  approach  to 
adaptive  detection  in  nonstationary  noise  and  clutter  background  is  to  compare  the  processed 
target  signal  to  an  adaptive  threshold.  In  the  cell-averaging  CFAR  processing,  an  estimate 
of  the  background  noise  from  the  leading  and  the  lagging  reference  windows  is  used  to  set 
the  adaptive  threshold.  A  threshold  multiplier  (or  scaling  factor)  is  used  to  scale  the 
threshold  to  achieve  the  desired  probability  of  false-alarm. 

In  the  first  part  of  this  report,  we  have  proposed  two  modified  cell-averaging  detectors  for 
multiple  target  situations.  The  first  one  is  a  weighted  cell-averaging  CFAR  detector. 


20  DISTRIBUTION /Avail  ABU  ITY  OF  ABSTRACT 

□  unclassified/unlimited  B  same  as  rpt 

22a  NAME  OF  RESPONSIBLE  INDIVIDUAL 

Vincent  C.  Vannicola 

DO  Form  1473.  JUN  86 

ATP  FOP'-E,  Sfl  1V29-4-88-  138  r  Wt) 


21  ABSTRACT  SECURITY  CLASSIFICATION 

□  dtic  USERS  UNCLASSIFIED 

22b  TELEPHONE  (Include  Area  Code)  22c  OFFICE  SYMBOL 

(315)  330-4437  RADC  (OCTS) 

Previous  editions  are  obsolete  SECURITY  CLASSIFICATION  OF  THIS  PAGE 

UNCLASSIFIED 


.nV.Iv  „<v 


rasa m 


m 


l&p 


UNCLASSIFIED 


^  WCA-CFAR,  where  weighted  leading  and  lagging  reference  windows  are  used  to  obtain  the 

adaptive  threshold.  The  second  is  a  cell-censored,  cell  averaging  CFAR  processor  where 
a  predetermined  fixed  threshold  is  used  to  eliminate  those  cells  that  may  contain 
interference.  _ _ .  >'  f.  i 

In  the  second  part  of  the  report,  the  theory  of  distributed  CFAR  detection  with  data 
fusion  is  developed.  First,  a  system  consisting  of  n  CA-CFAR  detectors  with  data 
fusion  is  considered.  The  overall  system  is  optimized  so  that  the  overall  probability 
of  detection  is  maximum  while  the  overall  probability  of  false-alarui  xs  iixed  at  uie 
desired  value.  Next,  CFAR  detection  with  multiple  background  estimators  and  a  data 
fusion  center  is  studied.  Finally,  adaptive  CFAR  detection  with  multiple  detectors 
for  different  network  topologies  is  considered.  Two  topologies,  namely,  a  parallel 
and  a  tandem  topology  are  investigated.  The  overall  systems  are  optimized  so  that 
the  probability  of  detection  is  maximum  while  CFAR  is  achieved. 


Accession  For  . 

8TIS  GRA&I 

c/ 

DTIC  TAB 

□ 

Unannounced 

□ 

1  . 

- - — 

Distribution/ 

Availability 

Codes 

Dlst 


lAvall  and/or 
Special 


\ 


UNCLASSIFIED 


TABLE  OF  CONTENTS 


LIST  OF  ILLUSTRATIONS . . . - . 

LIST  OF  TABLES . - . . . - 

CHAPTER  I.  INTRODUCTION . 

1.1  Constant  False-Alarm  Rate 

(CFAR)  Processors . - . 

1.1- 1  Background  and 

Literature  Review . 

1.1- 2  System  Description, 

Notation,  and  Terminology- 

1.2  Detection  with  Distributed 

Sensors  and  Data  Fusion-- . 

1.3.  Dissertation  Organization . 

CHAPTER  II.*  A  WEIGHTED  ADAPTIVE  CELL- AVERAGING 
CFAR  DETECTOR  FOR  MULTIPLE  TARGET 
SITUATIONS . 

2.1  Introduction . 

2.2  Preliminaries . . - . 

2.3  Cell-Averaging  CFAR  Processors - 

2.3- 1  CA-CFAR  Detector . 

2.3- 2  GO-CFAR  Detector-- . 

2.3- 3  SO-CFAR  Detector . 

2.3- 4  Threshold  Multipliers . 

2.4  Weighted  Cell-Averaging  CFAR 

Detector . . 

2.5  Summary  and  Conclusions . . 

CHAPTER  III.  A  CELL- CENSORED  CFAR  DETECTOR  FOR 

MULTIPLE  TARGET  SITUATIONS . 


3.1  Introduction 


CHAPTER  IV. 


Cell-Censored  CFAR  Detection- 


3.3  Summary  and  Conclusions- 


CA-CFAR  DETECTION  WITH  DISTRIBUTED 
RADARS  AND  DATA  FUSION . 


Introduction- 


CHAPTER 


CHAPTER 


Distributed  CFAR  Detection  with 
Data  Fusion . 


4.2-1  ’’AND"  Fusion  Rule- 


4.2-2  "OR"  Fusion  Rule 


4.2- 3  Example . - . 

4.2- 4  A  Comment  on  the  Perform¬ 

ance  of  the  "AND"  and  "OR" 
Fusion  Rules . 

4.3  Optimization  of  the  Overall 

System . . 

4.4  Summary  and  Discussion . 

ADAPTIVE  CFAR  DETECTION  WITH  MULTIPLE 
BACKGROUND  ESTIMATORS . 


Introduction . 

Adaptive  CFAR  Detection  with 
Multiple  Background  Estimators - 

5.2-1  CA-CFAR  Detection . 


GO-CFAR  Detection- 
SO-CFAR  Detection- 


Example-- . - . 

An  Extension . - . 

Summary  and  Discussion- 


ADAPTIVE  CFAR  DETECTION  WITH  DIFFERENT 
NETWORK  TOPOLOGIES . 


'vVw.'V 


6 . 1  Introduction- 


6.2  System  Optimization  and 

Performance  Analysis . 

6.2-1  Adaptive  CFAR  Detection 
with  Parallel  Network 
Topology . 


\>Nv- 

m  *  •  ' 

1  o 

'Wfw 


J  Vi 


Page 


6.2-2  Adaptive  CFAR  Detection 
with  Tandem  Network 
Topology .  127 

6.3  Example .  131 

6.4  A  Comparison .  136 

6.5  Summary  and  Conclusions .  138 

CHAPTER  VII.  SUMMARY  AND  SUGGESTIONS  FOR  FUTURE 

RESEARCH  WORK .  140 

7.1  Summary .  140 

7.2  Suggestions  for  Future  Research 

Work . - .  143 

APPENDIX  A .  144 

REFERENCES .  152 


* 

h 

h 

V 


* 


* 

t 

t 

< 

I 

i 

< 


111 


I  iT*  VWY  WUVWW lTVWWWVS 


Figure  No, 


LIST  OF  ILLUSTRATIONS 


Title 


Page 


Effect  of  the  Noise  Power  Increase 
on  the  Probability  of  False-Alarm 
for  a  Fixed  Threshold. 

Design  Pp  =  10  ® - - - 

Cell-Averaging  CFAR  Detector . . 

Variation  of  Noise  Power  as  a 
Function  of  Range . - . . 


Distributed  Sensor  System  with 
Central  Computation . . 


Distributed  Sensor  System  with 
Data  Fusion . . 


Weighted  Cell -Averaging  CFAR 
Detector . 


Detection  Probability  for  CA,  GO,  SO 
and  WCA-CFAR  Detectors  for 

N-4,  PF=10'4  and  I=S-  10- . - . 

Detection  Probability  for  CA,  SO,  GO 
and  WCA-CFAR  Detectors  for  N=4, 

P  =10'4  and  I  =  S . ---■ 

F 


Detection  Probability  CA,  GO,  SO 
and  WCA-CFAR  Detectors  for  N=4, 

P  =10'4  and  I  =  S  x  10 . . 

r 

Detection  Probability  for  WCA-CFAR 

-  4 

Detector  when  N=4,  P  =10  , 

I=S  10,  I=S ,  and  I*S  x  10 . 

Detection  Probability  for  WCA-CFAR 

Detector  when  Pp=10  4,  I=S  •=■  10, 
N=4  and  N=8 . . . 


m 


wf 


m 


s/v> 


.  ►  * * ,« * *  -*»• . » **  n" 


Figure  No . 


1  li* J*1  V J# Jk’ Jl* Jt* .1  iMi*  >»■.+»» 

Km 


._  Title 

Detection  Probability  for  WCA-CFAR 

Detector  when  P  =10  4  I=S ,  N=  4 
and  N=8 . ? . 


Page 


Detection  Probability  for  WCA-CFAR 

-  4 

Detector  when  PR=10  ,  I=S  x  10, 

N=4  and  N=8 . . . 


Cell-Censored  CFAR  Detector- 


Performance  of  the  Cell -Censored  CFAR  Detector 
PF=10‘5,  PQ=0 . 5  and  I=S  t  10 . . 

A  Performance  Comparison  of  the  Cell- 
Censored  (CC)  CFAR  Detector  and  the  CA-CFAR 

Detector.  P  =10"5,  PR  =  10'6,  Pn=0.5 

and  I=S  t  10?-- . ?!--- . - . 

A  Performance  Comparison  of  the  Cell - Censored 
(CC)  CFAR  Detector  and  the  CA-CFAR  Detector. 

P  =10'6,  pj  =10"3,  PA=0 . 5  and  I=S  t  10 . 

F  ^1  O 

A  Performance  Comparison  of  the  Cell- Censored 
(CC)  CFAR  Detector  and  the  CA-CFAR  Detector. 

Pp=10'5,  Pp  =10* 1 ,  PQ=0 . 5  and  I-S  t  10 . 

Distributed  CA-CFAR  Detection  with  Data 
Fusion . 

Distributed  CA-CFAR  Detection  and  Data 
Fusion  with  Two  Detectors . . 


Performance  Comparison  of  the  Two  Sensor 
System  (DS)  Using  the  "AND"  Fusion  Rule 
with  a  Single  Sensor  System  (SS) 

Pp=10'4,  N1  =  6,  and  N2  =  4 . . 

Performance  Comparison  of  the  Two  Sensor 
System  (DS)  Using  the  "OR"  Fusion  Rule 
with  a  Single  Sensor  System  (SS) . 

Pp=10"4,  N1  =  6  and  N2  =  4 . - . . 

Performance  of  the  Two  System  System  for 
the  "AND"  and  the  "OR"  fusion  Rule. 

Pp=10  4,  N^  =  6  and  N2  =  4 . - . . 


VFF, 

O  kXV. 


vV  V 
# 


‘V'V1 


s 


Figure  No. 


Pp  Versus  Pp  for  the  Overall 


Optimum  System.  S=5  dB- 


Pp  Versus  Pp  for  the  Overall 


Optimum  System  S-^IO  dB . 

CFAR  Processor  with  Multiple 
Background  Estimators . . 


Cell-Averaging  CFAR  Processor  with 
Multiple  Background  Estimators . . 

A  Two-Sensor  CFAR  Processor  with  One 
Background  Estimator . . 


Performance  of  the  CA-CFAR  Detectors 

for  P  =10'4  and  N=4 . . 

r 

Performance  of  the  GO-CFAR  Detectors 

for  Pc=10'4  and  N=4 - 

r 

Performance  of  the  CA,  GO  and  SO  CFAR 

Detectors  with  a  Background  Estimator 

for  P  =10-4  and  n=  4 . 

F 

Hybrid  CFAR  Processor  with  Data 
Fusion . 

Hybrid  System  with  Two  Background 
Estimators  and  a  Data  Fusion  Center- - 

Performance  of  the  Hybrid  System 
with  Two  Background  Estimators  (HS) 
and  the  Two-Sensor  Svstem  of 
Figure  4.2  (DS) ,  Pp=10-*,  N,-6 
and  N2=4 . - . - . 

Parallel  Network  Topology  for 
Adaptive  CFAR  Detection - - 


Tandem  Network  Topology  for 
Adaptive  CFAR  Detection - 


A  Two  Sensor  Network  Topology  for 
Adaptive  CFAR  Detection . . 


Different  Network  Topologies  for 
Two-Sensor  CFAR  Detection  Systems - 


.'y 

v v  v  v  v  v  v  v  v  V 


LIST  OF  TABLES 


Table  No 


Title 


’AND"  Fusion  Rule- 
’OR"  Fusion  Rule- 


Performance  of  Adaptive  CFAR 
Detection  with  Two-Sensor  System 
and  a  Single  Sensor  System. 

Pp  =  10'5,  N=6 ,  Nx  =  4 . 

CA-CFAR  Performance  of  the  Two- 
Sensor  Systems  of  Figure  6-4 . 


g 


CHAPTER  I 


INTRODUCTION 


1.1  Constant  False  Alarm  -  Rate  fCFAR)  Processors 


1.1-1  Background  and  Literature  Review 

In  practical  radar  signal  detection  systems,  the 
problem  is  to  automatically  detect  a  target  in  a  nonstation¬ 
ary  noise  and  clutter  background  while  maintaining  a  constant 
probability  of  false-alarm.  "Clutter"  is  the  term  applied 
to  any  unwanted  radar  signal  from  scatterers  that  are  not  of 
interest  to  the  radar  user.  Examples  of  unwanted  echoes,  or 
clutter,  in  radar  signal  detection  are  reflections  from 
terrain,  sea,  rain,  birds,  insects,  chaff,  etc.  Classical 
detection  using  a  matched  filter  receiver  and  a  fixed 
threshold  is  not  applicable  due  to  the  nonstationary  nature 
of  the  background  noise.  In  fact,  a  small  increase  in  the 
total  noise  power  results  in  a  corresponding  increase  of 
several  orders  of  magnitude  in  the  probability  of  false 
alarm.  As  shown  in  Figure  1-1,  for  a  design  probability  of 
false-alarm  of  10  ,  an  increase  of  only  3  dB  in  signal-to- 

noise  ratio  causes  the  actual  fclse-alarm  probability  to 

4 

increase  by  almost  10  .  This  undesirable  increase  in  the 
number  of  false  alarms  causes  the  data  processing  equipment, 
either  a  human  operator  or  a  digital  computer,  to  saturate. 


SMtwS* 


Therefore,  adaptive  threshold  techniques  are  needed  to  main¬ 


tain  a  constant  false -alarm  rate  (CFAR) . 


One  approach  to  adaptive  detection  in  nonstationary 


noise  and  clutter  background  is  to  compare  the  processed 


target  signal  from  the  test  cell  to  an  adaptive  threshold 


obtained  from  the  mean  level  of  clutter  plus  noise  over  ad¬ 


jacent  range  and/or  Doppler  Cells  [1-3].  The  conventional 


cell-averaging  constant  false-alarm  rate  detector  is  shown 


in  Figure  1-2.  In  the  cell -averaging  constant  false-alarm 


rate  detector,  CA-CFAR,  the  adaptive  threshold  is  obtained 


from  the  arithmetic  mean  of  the  reference  cells  [4].  Finn  [2] 


showed  that  in  a  homogeneous  Gaussian  noise  background,  where 


the  noise  samples  obtained  from  the  range  cells  are  identically 
distributed,  the  CA-CFAR  detector  performs  well  and  its  per¬ 


formance  approaches  that  of  the  Neyman- Pearson  detector  as 


the  number  of  range  cells  is  increased  to  infinity.  In 


nonhomogeneous  background,  which  may  be  caused  by  clutter 


edges  and  chaff,  the  adaptive  threshold  setting  is  seriously 


affected  resulting  in  a  degradation  of  the  performance.  In 


the  situation  where  the  transition  from  a  clear  to  a  clutter 


environment  is  not  relatively  smooth,  it  is  assumed  that 


the  total  noise  power  density  as  a  function  of  range  can  be 


represented  by  the  step  function  as  shown  in  Figure  1-3, 


Two  cases  may  be  encountered  in  this  severe  clutter  environ¬ 


ment  [1,2].  In  the  first  case,  the  cell  under  test  is  in 


the  clear  region  but  a  group  of  the  reference  cells  are 


immersed  in  the  clutter.  This  results  in  a  higher  adaptive 


»w*  r*  i**,  n,  n. 


threshold  and  the  probabilities  of  detection  and  false-alarm 
are  reduced.  This  is  also  known  as  the  masking  effect.  In 
the  second  case,  if  the  cell  under  test  is  immersed  in  the 
clutter  but  some  of  the  reference  cells  are  in  the  clear 
region,  the  threshold  is  relatively  low  and  the  probability 
of  false-alarm  is  increased  [1,2, 5, 6].  Jamming  and  the 
clutter  edges  mentioned  above  significantly  affect  the  power 
level  as  a  function  of  range  which  results  in  an  intolerable 
increase  in  the  probability  of  false-alarm.  This  situation 
is  frequently  encountered  in  search  radar  design.  To  control 
this  intolerable  increase  in  the  probability  of  false-alarm, 
Hansen  and  Sawyers  [5,6]  proposed  the  greatest-of-selection 
logic  in  cell-averaging  constant  false-alarm  rate  detector, 
GO-CFAR.  A  comparison  between  the  two  CFAR  detectors,  CA-CFAR 
and  GO-CFAR,  using  square  law  detection  of  Swerling  I  targets 
was  treated  by  Moore  and  Lawrence  [7].  The  clutter  map  CFAR 
approach  can  also  be  used  for  nonhomogeneous  background.  In 
this  technique  proposed  by  O’Donnell,  Muehe,  and  Labitt  [8], 
the  detector  output  of  each  resolution  cell  is  averaged  over 
several  scans.  An  analysis  of  the  clutter  map  CFAR  technique 
has  been  performed  by  Nitzberg  [9],  to  obtain  the  probability 
of  detection. 

Trunk  [10],  while  studying  the  target  range  resolu¬ 
tion  of  some  adaptive  threshold  detectors,  showed  that  targets 
cannot  be  resolved  by  a  CA-CFAR  detector  if  another  target  lies 
within  the  reference  cells  of  the  other.  He  proposed  the  smal- 
lest-of-selection  logic  in  cel 1  -  averaging  constant  false-alarm 


rate  detector,  SO-CFAR.  The  SO-CFAR  detector  is  less  sensi¬ 
tive  to  the  detection  loss  than  the  CA-CFAR  unless  the  number 
of  resolution  cells  is  relatively  large.  In  other  words,  for 
a  small  number  of  reference  cells,  the  detection  loss  of  the 
SO-CFAR  is  relatively  large  as  compared  to  the  ideal  detector. 
Closely  separated  targets  may  be  encountered  in  dense  civilian 
or  military  target  environments.  In  such  situations,  capture 
effect  may  occur  where  the  probability  of  detection  rises 
slowly  and  reaches  one  asymptotically.  Finn  and  Johnson  [1] , 
and  Rickard  and  Dillard  [11]  studied  this  problem  of  two 
closely  separated  targets  and  showed  the  existence  of  the 
capture  effect  for  all  single-pulse  Swerling  target  models 
(cases  1,2,  and  4).  The  SO-CFAR  detector  proposed  in  [10] 
eliminates  the  capture  effect  that  exists  in  the  CA-CFAR  and 
GO-CFAR  detectors. 

Additional  work  on  CFAR  processors  in  situations 
where  one  or  more  targets  lie  in  the  set  of  reference  cells, 
i.e.,  in  multiple  target  situations,  is  reported  in  the 
literature  [12-16].  In  this  situation,  the  interfering 
targets  raise  the  threshold,  thereby  drastically  reducing  the 
probability  of  detection.  The  noise  samples  of  the  reference 
cells  are  not  identically  distributed  in  this  case.  There¬ 
fore,  a  mathematical  representation  of  the  system  model 
is  more  complicated  which  makes  the  performance  analy¬ 
sis  of  the  systems  more  involved.  McLane  et  al.  [12]  proposed 
a  threshold  control  technique  for  the  CA-CFAR  detector  based 
on  a  priori  information  about  location  that  could  be  supplied 


by  the  radar's  tracking  system.  A1  Hussaini  and  Ibrahim  [13] 
presented  a  generalization  and  extension  of  the  technique  in 
[12]  for  the  CA-CFAR,  GO-CFAR,  and  SO-CFAR  detectors.  When 
the  interfering  target  is  in  one  of  the  reference  cells, 

Weiss  [14]  showed  that  the  detection  for  a  GO-CFAR  detector 
is  extremely  poor.  To  alleviate  this  problem,  he  suggested 
the  use  of  the  SO-CFAR  technique.  Rohling  [15]  introduced  an 
order  statistic  based  estimation  technique  to  obtain  CFAR  in 
multiple  target  situations.  Finally,  Ritcey  [16]  presented 
the  performance  analysis  of  the  censored  mean-level  detector 
(CMLD) ,  which  is  an  alternative  to  the  mean-level  detector 
(MLD)  considered  previously  [1-15].  He  obtained  the  expres¬ 
sions  for  the  probability  of  detection  of  the  C MLD  in  a 
multiple- target  environment  when  a  fixed  number  of  Swerling  II 
targets  are  present.  The  CMLD  achieves  robust  detection  per¬ 
formance  in  a  multiple- target  environment  by  censoring 
several  of  the  largest  samples  of  the  maximum  likelihood 
estimate  of  the  background  noise  level. 

1.1-2.  System  Description,  Notation,  and 
Terminology 

In  this  subsection,  we  describe  the  operation  of  the 
cell-averaging  CFAR  detector,  we  define  some  notation,  and 
present  some  mathematical  background  that  will  be  used 
throughout  the  report. 


The  conventional  cell-averaging  CFAR  processor  is 
shown  in  Figure  1-2.  As  shown,  the  output  from  the  square 


law  detector  is  fed  into  a  tapped  delay  line  torming  the 
reference  cells.  The  delay,  x,  between  the  taps  is  approxi¬ 
mately  equal  to  the  transmitted  pulsewidth.  To  avoid  any 
signal  energy  spill  from  the  test  cell  into  the  directly 
adjacent  range  cells,  which  may  affect  the  clutter  power 
estimation,  the  adjacent  cells,  called  guard  cells  (usually 
one  or  two  on  each  side  of  the  cell  under  test),  are  com¬ 
pletely  ignored.  The  statistics  of  the  reference  windows  U 

N 

and  V  are  obtained  from  the  sum  of  the  y  leading  cells  and 
the  ^  lagging  reference  cells  respectively.  Thus,  a  total 
of  N  noise  samples  are  used  to  estimate  the  background  en¬ 
vironment.  The  noise  samples  are  assumed  to  be  statistically 
independent.  The  reference  windows  U  and  V  are  combined 
according  to  some  selection  logic  to  obtain  the  estimate  of 
the  clutter  power  level  Q.  To  maintain  the  probability  of 
false-alarm,  Pp,  at  the  desired  value,  the  adaptive  thresh¬ 
old  is  multiplied  by  a  scaling  factor  called  the  threshold 
multiplier  T^.  It  should  be  noted  that  for  the  same  value 
of  Pp,  the  threshold  multiplier  is  not  the  same  for  different 
selection  logics.  The  product  T^Q  is  the  resulting  adaptive 
threshold.  The  output,  qQ,  from  the  center  tap  is  then  com¬ 
pared  with  the  threshold  in  order  to  make  a  decision. 

We  assume  that  the  target  at  the  test  cell  (center 
tap),  called  the  primary  target,  is  a  slowly  fluctuating 
target  model  of  Swerling  type  I.  The  signal-to-noise  ratio, 
SNR,  of  the  primary  target  is  denoted  by  S.  We  further 
assume  that  the  total  background  noise  is  white  Gaussian. 
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Since  both  the  noise  and  Rayleigh  targets  have  Gaussian 

quadrature  components,  the  output  of  the  square-law  detector 

has  an  exponential  probability  density  function  [17].  If 

2 

the  noise  variance  is  a  ,  then  the  conditional  probability 
density  function  of  the  output  of  the  detection  cell,  qQ,  is 
given  by 


PQ  ’ 

o  1  X 


-  qQ/ 2cr  (1  +  S) 


2a  (1+S) 


1  'V2a' 


,  for  hypothesis 
,  for  hypothesis  H 


(1-D 


The  hypothesis  Hq  represents  the  case  of  noise  alone,  while 
hypothesis  represents  the  noise  plus  target  signal  case. 
The  probability  of  detection,  P^,  where  N  may  represent  CA, 
GO,  or  SO,  is  given  by 
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PD  *  |  PrfQo  >  TNq  |  Q,  Hj)  I>|*(q)  dq  (1-2) 
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Pq(q)  denotes  the  probability  density  function  of  the  adap¬ 


tive  threshold  and 


oo 

prCQ0>TNqIQ,Hi)  =  pq  |  h.  ^qo  lHl^dqo  =  expt"  T~  7T] 


(1-3) 


We  easily  see  that  Pp  is  directly  obtained  from  P^  by  just 

setting  the  target  SNR,  S,  to  zero,  i.e. 

00 

Pp  -  )  Pr(Q0  >  TNq  |  Q,  Hq)  pj*(q)  dq 


'TNq/2a  N 
e  N  •  Pq (q)  dq 


d-4) 


.  •  .  X 


For  simplicity,  we  normalize  the  noise  power  and  the  output 


from  the  test  cell  and  divide  by  2 a  .  Then,  equations  (1-1) 


and  (1-3)  become 


1  -q0/d+s) 

yyg  e  ,  for  hypothesis  H, 


(1-5) 


,  for  hypothesis  H 


pr(Q0  >  TNq|Q,  =  e 


(1-6) 


When  there  is  no  interference  in  the  reference  cells,  then 
all  the  noise  samples  are  identically  distributed.  Since  these 
noise  samples  are  statistically  independent,  the  reference 
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windows  U  and  V  have  a  chi-square,  x  >  distribution  of  2 (y) 

degrees  of  freedom  (see  Appendix  A  for  details),  i.e., 
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The  corresponding  cumulative  distributions  are 


Fy(u)  =  1  -  e 
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Fv(v)  =  1  *  e  l  -t  vJ  ,  v  >_  0 
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(1-10) 


The  above  results  will  be  used  throughout  this  dissertation. 


Next,  we  discuss  distributed  detection  and  data  fusion. 


1.2.  Detection  with  Distributed  Sensors  and  Data  Fusion 

The  use  of  multiple  sensors  with  data  fusion  is  widely 
increasing  in  surveillance  systems.  One  of  the  main  goals 
of  using  multiple  sensors  is  to  improve  system  performance 
such  as  reliability  and  speed.  For  surveillance  systems 
requiring  a  large  area  of  coverage  and/or  a  large  number  of 
targets  under  consideration,  multiple  sensors  are  used.  In 
such  systems,  complete  observations  can  be  transmitted  by  the  . 
sensors  to  a  central  processor  for  data  processing  as  shown 
in  Figure  1-4.  This  requires  a  large  communication  bandwidth 
which  may  not  be  available.  Thus,  due  to  the  constraints  on 
the  bandwidths  of  the  communication  channels,  distributed 
signal  processing  with  a  data  fusion  center  is  preferred  in 
many  situations.  In  such  distributed  detection  systems,  some 
processing  of  the  signal  is  done  at  each  sensor  which  then 
sends  partial  results  to  the  data  fusion  center,  as  shown 
in  Figure  1-5.  These  partial  results  are  combined  according 
to  a  suitable  data  fusion  rule  to  yield  the  desired  global 
result.  In  our  case,  the  partial  results  are  the  decisions 
from  the  individual  detectors,  D^,  i=l,2,...,n,  where 

e  {0,1}.  The  D^'s  are  combined  to  yield  a  final  decision, 
Dq,  which  may  again  be  a  zero  or  a  one. 

Some  work  on  distributed  detection  has  been  reported 
in  the  literature  [18-31].  Tenney  and  Sandell  [18],  extended 
the  classical  Bayesian  theory  to  the  case  of  distributed 
sensors  for  the  binary  hypothesis  testing  problem.  Sadjadi 
[19]  generalized  the  results  of  [18]  to  the  case  of  multiple 


Phenomenon  H 


Data  Fusion  Center 


hypotheses.  Ekchian  and  Tenney  [20]  considered  various 
distributed  sensor  network  topologies  and  solved  some  associ¬ 
ated  Bayesian  hypothesis  testing  problems.  Teneketzis  [21, 
22]  has  solved  the  decentralized  quickest  detection  problem 
and  the  decentralized  Wald  problem.  Srinivasan  [23]  has 
presented  some  results  on  distributed  radar  detection. 

Kushner  and  Pacut  [24]  conducted  a  simulation  study  on  the 
effects  of  specific  prior  probabilities  and  parametric  de¬ 
pendencies  on  the  decision  rule.  Chao  and  Lee  [25]  presented 
a  distributed  detection  scheme  based  on  soft  local  decisions. 
Conte,  D 'Addio,  Farina  and  Longo  [26]  studied  the  design  and 
performance  evaluation  of  optimum  and  suboptimum  multistatic 
radar  receivers.  Chair  and  Varshney  [27]  derived  the  optimum 
fusion  rule  at  the  data  fusion  center  for  a  Bayesian  detec¬ 
tion  problem  with  distributed  sensors.  Hoballah  and  Varshney 
[28]  solved  the  problem  of  Neyman- Pearson  detection  with  dis¬ 
tributed  radars.  In  [29],  they  devised  an  information  theo¬ 
retic  formulation  to  the  distributed  detection  problem.  More 
work  on  distributed  detection  and  data  fusion  has  been  done 
by  Varshney  and  his  coworkers  [30,31]. 


1  •  3  .  Report  Organization 

Tn  this  report,  two  major  topics  are  considered. 
The  first  part  deals  with  CFAR  detection  in  multiple  target 
situations.  In  Chapter  two,  we  present  a  weighted  cell¬ 
averaging  CFAR  detector  for  multiple  target  situations.  The 
leading  reference  window  and  the  lagging  reference  window  are 
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weighted  and  added  to  yield  the  adaptive  threshold.  The 
relative  weights  assigned  to  the  reference  windows  depend 
upon  the  level  of  interference.  These  weights  are  selected 
such  that  the  desired  probability  of  false-alarm  is  maintained 
while  the  probability  of  detection  is  maximized.  The  per¬ 
formance  of  the  weighted  cell-averaging  CFAR  detector  is 
analyzed  and  compared  to  the  performance  of  the  conventional 
cell-averaging,  the  greatest-of- selection  cell -averaging  and 
the  smallest-of- selection  cell-averaging  CFAR  detectors. 

In  Chapter  three,  we  present  a  cell-censored  CFAR  de¬ 
tector.  Each  range  cell  is  first  compared  to  a  predetermined 
fixed  threshold,  where  a  decision  about  the  presence  or  absence 
of  an  interfering  target  is  made.  The  cells,  where  the  presence 
of  an  interfering  target  is  decided,  are  not  used  while  forming 
the  adaptive  threshold.  This  censoring  scheme  eliminates  the 
cells  with  interfering  targets  which  may  otherwise  raise  the 
threshold  and  therefore,  lower  the  probability  of  detection. 

The  performance  of  the  cell -censored  CFAR  detectr  is  also 
studied. 

The  second  part  of  this  report  deals  with  cell¬ 
averaging  CFAR  detection  using  multiple  sensors  and  data  fusion. 
In  Chapter  four,  we  study  a  system  consisting  of  n  CA-CFAR 
detectors  and  a  data  fusion  center.  Given  the  fusion  rule, 
first  we  obtain  the  optimum  threshold  multipliers  at  the  indi¬ 
vidual  CFAR  detectors  by  maximizing  the  overall  probability  of 
detection  at  the  data  fusion  center,  while  the  overall  proba¬ 
bility  of  false-alarm  is  maintained  at  the  desired  value. 
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Next,  we  consider  the  optimization  of  the  overall  system,  i.e., 
we  derive  both  the  optimum  fusion  rule  as  well  as  the  optimum 
threshold  multipliers  at  individual  CA-CFAR  detectors. 


In  Chapter  five,  we  consider  a  system  where  multiple 
background  estimators  are  used  for  the  estimation  of  the  back¬ 
ground  noise.  These  estimates  are  transmitted  to  the  CFAR 
detector.  This  CFAR  detector  computes  its  own  estimate  of  the 
background  noise  and  combines  it  with  the  received  estimates 
to  yield  the  adaptive  threshold.  The  performance  of  this 
system  is  studied  and  an  extension  is  also  proposed. 

In  Chapter  six,  we  consider  adaptive  CFAR  detection 
for  two  distributed  sensor  network  topologies.  Specifically, 
we  consider  the  parallel  and  the  tandem  network  topologies. 

In  this  chapter, the  compressed  data  transmitted  amongst  the 
detectors  is  assumed  to  be  in  the  form  of  decisions  instead 
of  estimates.  The  overall  systems  are  optimized  to  yield 
the  maximum  probability  of  detection  for  a  fixed  probability 
of  false-alarm.  The  performance  of  the  systems  is  also 
analyzed. 

Finally,  in  Chapter  seven,  we  present  a  summary  along 
with  a  discussion.  Some  suggestions  for  future  research  work 
are  also  presented. 
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CHAPTER  II 


A  WEIGHTED  ADAPTIVE  CELL-AVERAGING  CFAR 
DETECTOR  FOR  MULTIPLE  TARGET  SITUATIONS 


2.1  Introduction 


In  radar  signal  detection,  the  problem  is  to  auto¬ 
matically  detect  a  target  in  a  nonstationary  noise  and 
clutter  background  while  maintaining  a  constant  probability 
of  false-alarm.  Classical  detection  using  a  matched  filter 
receiver  and  a  fixed  threshold  is  not  applicable  due 
to  the  nonstationary  nature  of  the  background  noise.  A 
small  increase  in  the  total  noise  power  results  in  a  cor¬ 
responding  increase  of  several  orders  of  magnitude  in  the 
probability  of  false-alarm.  Therefore,  adaptive  threshold 
techniques  are  needed  to  maintain  CFAR. 

In  this  chapter,  we  propose  and  analyze  a' modified 
cell -averaging  detector  for  multiple  target  situations. 

The  leading  and  the  lagging  windows  are  weighted  in 
accordance  with  the  level  of  interference.  Optimum 
weights,  which  maintain  a  constant  probability  of  false- 
alarm  and  maximize  the  probability  of  detection, are  obtained 
For  the  Swerling  fluctuating  target  model  I,  we  obtain  the 
expression  for  the  probability  of  detection  for  our  weighted 
cell- averaging  CFAR  detector.  In  Section  2.2,  we  present 
some  background  material  and  assumptions.  In  Section  2.3, 
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we  derive  the  expression  for  the  probability  of  detection 
for  the  CA-CFAR,  GO-CFAR,  and  SO-CFAR  detectors.  In  Section 
2.4,  we  evaluate  the  performance  of  the  weighted  cell¬ 
averaging  CFAR  detector,  WCA-CFAR.  Numerical  results  show¬ 
ing  the  performance  of  the  WCA-CFAR  detector  are  also  pre¬ 
sented.  In  Section  2.5,  we  provide  a  summary  along  with  a 
discussion  of  the  results. 


2 . 2  Preliminaries 

The  conventional  cell-averaging  CFAR  processor,  shown 
in  Figure  1-2,  has  been  described  in  Section  1.1-2.  There, its 
performance  in  the  presence  ot  interference  was  not  discussed. 
In  this  chapter,  we  consider  cell-averaging  CFAR  processors, 
in  the  presence  of  interference,  i.e.,  in  multiple  target 
situations.  We  assume  that  the  target  at  the  test  cell 
(center  tap),  called  the  primary  target  is  a  slowly  fluctu¬ 
ating  target  model  of  Swerling  type  I.  The  signal-to-noise 
ratio,  SNR,  of  the  primary  target  is  denoted  by  S.  Without 
loss  of  generality,  we  assume  only  one  interfering  tar¬ 
get,  the  secondary  target,  in  one  of  the  taps  of  the  window  U. 
We  assume  that  this  interfering  target  is  also  fluctuating  in 
accordance  with  Swerling  target  model  I  with  an  SNR  I .  We 
further  assume  that  the  total  background  noise  is  white  Gaus¬ 
sian.  Since  the  noise  samples  of  the  reference  window  V  are 
identically  distributed,  the  probability  density  function  and 
the  cumulative  distribution  function  are  as  defined  in  (1-8) 
and  (1-10).  Due  to  interference,  the  noise  samples  of  U  are 
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(2-4) 


Using  Leibniz's  rule,  we  differentiate  equation  (2-4)  to 
obtain  the  probability  density  function  of  the  adaptive 
threshold  which  is  just  the  convolution  of  the  random  vari¬ 
ables  U  and  V.  That  is, 

00 

PqA(q)  =  pu(u)  *  pv(v)  =  j  PyCu)  py(q-u)  du  (2-5) 

-  OO 

where  *  denotes  convolution  and  Leibniz's  rule  is  defined 
as  follows.  If 


b  (u) 

G(u)  =  H(x,  u)  dx 

a  (u) 


(2-6. a) 


then, 
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(2-6. b) 

Using  the  fact  that  U,V  0  and  U+V  =  Q,  equation  (2-5)  becomes 

CA  fq 

pq  (q)  =  ^(u)  py(q-u)  du  .  (2-7) 


(2-7) 


Substituting  equations  (2-1)  and  (1-8)  into  equation  (2-7), 
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2.3-3  SO-CFAR  Detect or 


In  this  case,  the  adaptive  threshold,  Q,  is  the  miminum 
of  the  reference  windows  (J  and  V,  i.e., 


Q  =  MIN (U,V) 


(2-19 


The  distribution  function  of  Q  is 


Fq  (q)  =  Pr  (Min (U, V)  <  q} 

=  Pr  (U  <  q  or  V  <  q}  (2-20] 


Using  the  complementary  function,  we  can  write 


1  -  Fq  Cq)  =  Pr  {min(U,Vj  >  q} 

=  Pr  {U  >  q  and  V  >  q} 


=  Pr  {U  >  q }  •  Pr{V  >  q} 


=  [1  -  Fu(q)][l  -  Fv(q)] 

Therefore,  the  cumulative  distribution  of  the  adaptive 
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threshold  is 


fq (q)  »  Fu(q)  +  Fy(q)  -  Fu(q)Fv(q) 


(2-22) 


Taking  the  derivative  of  equation  (2-22)  with  respect  to  q, 
we  obtain  the  first  order  probability  density  function  of  the 
adaptive  threshold  Q,  i.e.. 


Pq0 (q )  =  Py(q)  +  Pv(q)  -  Py(q)Fv(q)  -  Fu(q)Py(q) 

-  Pu(q)  [1  -  Fv(q)]  +  Py(q)  [1  -  F^q)]  (2-23) 


Substituting  equations  (1-8),  (1-10),  (2-1),  and  (2-2)  into 
equation  (2-23),  the  density  of  the  adaptive  threshold  becomes 
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Now,  substituting  equations  (2-24)  and  (1-6)  into  equa¬ 
tion  (1-2),  the  probability  of  detection  for  the  SO-CFAR 


detector  is  obtained  from 
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This  results  in 
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2.3-4  Threshold  Multipliers 


The  threshold  multipliers  for  the  three  schemes,  T^, 


Trf),  and  Tj,...,  are  all  different  since  they  depend  on  the 


selection  logic.  To  solve  for  these  threshold  multipliers, 
we  obtain  the  expression  for  the  probability  of  false-alarm 
of  each  detector.  Then,  we  set  the  probability  of  false- 


alarm  to  the  desired  value  v,  and  obtain  tca» tGq  and  Tgo’ 
from  equations  (2-27),  (2-28)  and  (2-29),  respectively,  by 
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2.4  Weighted  Cell- Averaging  CFAR  Detector 


The  proposed  weighted  cell-averaging  constant  false- 
alarm  rate  detector,  WCA-CFAR,  for  multiple  target  situ¬ 
ations  is  shown  in  Figure  2-1.  The  output  from  the  square 
law  detector  is  fed  into  a  tapped  delay  line  forming  the 
reference  cells.  The  set  of  the  leading  reference  cells 
form  the  reference  window  U,  while  the  set  of  the  lagging 
reference  cells  form  the  reference  window  V.  The  reference 
window  U  is  multiplied  by  the  weighting  coefficient  a,  while 
the  lagging  reference  window  V  is  multiplied  by  another 
weighting  coefficient,  g.  These  weighted  means  of  the 
leading  cells  and  the  lagging  cells  are  added  to  yield  the 
adaptive  threshold,  Q.  The  output,  qQ,  of  the  test  cell 
from  the  center  tap  is  compared  with  the  adaptive  threshold 
where  a  detection  decision  is  made. 

An  intuitive  justification  for  using  a  weighted  cell¬ 
averaging  approach  is  the  following.  While  computing  an 
estimate  of  the  background  noise  from  the  reference  cells, 
we  should  assign  more  weight  to  those  cells  that  do  not  con¬ 
tain  interference.  Thus,  if  interference  is  present  in  the 
leading  cells,  then  the  weight  a  assigned  to  the  reference 
window  U  should  be  smaller  than  the  weight  g  assigned  to  the 
reference  window  V.  The  relative  values  of  a  and  g  depend 
upon  the  level  of  interference  in  the  reference  cells.  It 
should  be  mentioned  that  a  and  g  should  be  such  that  the 
desired  probability  of  false-alarm  is  maintained. 
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Fig.  2-1.  Weighted  Cell- Averaging  CFAR  Detector. 


.  “  W  a"  W  *  «,  *  „  *  *  _  W  »  *  a  *  a  %  «  "  „  •  m'*m*  .  <,  «,  *  „  *  «  •  .  •  .  ■  »  *•  ■  '  »  ' ■  «  •  *  *'•  *  ■  «  *  ■  »  »  *■  a  *  »  *  .  »  k  *  »*  *•  *  «  ^  V  ■  ■*  *  *  *  -  *  ■  1 

_  trJ  IT  „  /.  r  .  /.  »"  »’  •_  if.  rnm  f  .  -  ,  /,  r.  r.  «,  /.  •»_  •  **,  »_  <  .  «  .  *  .  «  ,  ,  ■»,  *V  *  _  «_  *  «  ■r.  *  .  *  ,  *  I  '  ,  *  -  ■  *  1 

wV\VAV.Sw%  .A  aA -V^V/k^AV.%/' >  A  .  A ^  -VV«  A  ^  •  .<Vy^Vr>\V^\vV> 


re$i 

n 

rS8& 


M«jS» 

EwwS* 

rvv.N' 

lm: 


H 


mm 


We  assume  that  the  primary  target  at  the  test  cell 
has  a  signal- to-noise  ratio,  S,  and  the  secondary  target, 
or  interference  has  an  SNR,  I.  Without  loss  of  generality, 
we  assume  that  the  interference  is  in  the  leading  cells 
only.  The  primary  target  and  the  interfering  target  are 
assumed  to  be  fluctuating  in  accordance  with  the  Swerling 
target  model  I.  The  total  noise  power  in  the  background 
is  assumed  to  be  white  Gaussian.  The  reference  windows  U 
and  V  are  statistically  independent  random  variables.  Re¬ 
scaling  U  and  V  by  the  constants  a  and  g,  respectively,  we 
obtain  two  new  independent  random  variables  X  and  Y,  where, 


pY  (x)  =  -r-r  Pn(~)  ,  “  *  0  (2-30) 

X  |  ot  |  u  a 

and 

PyCy)  =  fjf  Pv(|)  ,  6*0  (2-31) 


Substituting  equations  (1-8)  and  (2-1)  into  equations 
(2-30)  and  (2-31),  we  obtain  the  probability  density  func¬ 
tions  of  X  and  Y  to  be 
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(2-32) 

(2-33) 


31 


The  adaptive  threshold,  Q,  of  the  weighted  cell -averaging 
detector  is  the  sum  of  X  and  Y,  i.e., 


Q  =  X  +  Y 


(2-34) 


Since  the  adaptive  threshold  Q  is  the  sum  of  the  two  inde¬ 
pendent  random  variables  X  and  Y,  the  probability  density 
function,  Pq(q)  is  the  convolution  of  X  and  Y,  i.e., 

00 

Pq(q)  =  px(x)  *  pY(y)  =  j  Vy)  px^-y)  dy  (2-35) 

-  00 

Using  the  constraints  (X,Y,  _>  0),  the  probability  density 
function  of  the  adaptive  threshold  becomes 


pq(q)  =  f  py(y)  px(q-y)dy  (2-36) 

0 

The  probability  of  detection  for  this  WCA-CFAR  detector 
is  defined  as 


PD  =  j  PrC%  >  IQ*  PQ^^  d(i 


(2-37) 


where 
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Pr(Qo>q|Q,  =  j  p  |  H_  (%'Hl)  d%  =  exp  ['  US] 
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(2-38) 


Now,  substituting  equations  (2-36)  and  (2-38)  into  (2-37), 
Pp  becomes 
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This  results  in 
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Using  the  identity  (2-9),  the  probability  of  detection  can 
be  rewritten  as 


PW  =  V 
rn  Y 
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Solving  the  double  integrals,  the  probability  of  detection 
is  given  by 


where  C,  C' ,  and  are  as  defined  in  (2-11). 

The  goal  of  our  WCA-CFAR  detector  is  to  obtain  the 

w 

weights  a  and  g  so  as  to  obtain  the  best  PD  while  keeping 

w 

the  false-aiarm  probability,  Pp,  to  be  a  constant  v.  This 
can  be  achieved  by  using  the  calculus  of  extrema,  for  which 


we  form  the  following  objective  function. 

J(«,g)  =  Pp(a,g)  +  ^  [Pp  Cot ,  B)  -  v] 


(2-42) 


where  v  is  the  desired  false-alarm  probability  and  £  is  the 
Lagrange  multiplier.  To  obtain  the  optimum  values  of  a  and 
g,  we  maximize  the  objective  function  J(a,g),  i.e.,  we  take 
the  derivatives  of  J(ot,g)  with  respect  to  a  and  g,  set  them 
equal  to  zero,  and  solve  for  a  and  g.  The  derivatives  are 
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We  obtain  £  from  equation  (2-43),  substitute  for  it  in 
equation  (2-44),  and  solve  for  6  in  terms  of  a.  We  obtain 


(1+S)  h2  -  h2 


where 
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Knowing  the  fact  that  the  probability  of  false-alarm  is 


(N/2) - 2 
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(2-46) 


we  substitute  for  3  from  equation  (2-45)  into  equation 
(2-46),  and  obtain  a  nonlinear  equation  in  terms  of  a  only, 
Using  an  iterative  procedure,  we  solve  for  a  and  B.  Since 
the  equations  are  nonlinear,  more  than  one  solution  may 
be  obtained.  Only  those  solutions  which  satisfy  the 
given  constraints  of  the  problem  are  kept. 


We  obtain  some  numerical  results  to  illustrate  the  perfor- 

W 

mance  of  the  WCA-CFAR  detector.  The  desired  Pp  is  assumed 
to  be  10 In  Figures  2-2,  2-3,  2-4,  we  compare  the  per¬ 
formance  of  the  WCA-CFAR  with  the  CA-CFAR,  GO-CFAR,  SO-LFAk 
detectors  in  terms  of  versus  SNR  curves  for  different 
interference  levels.  We  notice  that  the  WCA-CFAR  detector 
performs  better  than  the  other  three.  In  Figure  2-S,  we  show 
the  performance  of  the  WCA-CFAR  detector  as  a  function  of  the 
interference  level.  In  Figures  2-6,  2-7,  and  2-8  we  also  ob- 
interference  level.  As  expected,  in  Figures  2-6,  2-7,  and 
2-8  we  observe  that  the  performance  of  the  WCA-CFAR  improves 
when  the  number  of  reference  cells  is  increased. 


2.5  Summary  and  Conclusions 


In  this  chapter,  we  proposed  a  weighted  cell -averaging 
CFAR  detector  for  multiple  target  situations.  In  this  de¬ 
tector,  the  weighted  leading  and  lagging  reference  windows 
are  added  to  obtain  the  adaptive  threshold.  The  weights  are 
selected  so  that  CFAR  is  achieved  while  maximizing  the  proba¬ 
bility  of  detection.  It  results  in  a  smaller  weight  being 
assigned  to  the  window  where  interference  occurs  and  more 
weight  is  assigned  to  the  window  without  interference.  This 
weighting  procedure  prevents  the  adaptive  threshold  from  be¬ 
coming  too  high  and  resulting  in  a  severe  degradation  of  the 
probability  of  detection.  For  Swerling  target  model  I  embedded 
in  a  white  Gaussian  noise  of  unknown  level,  we  showed  that 
the  WCA-CFAR  detector  performs  better  than  the  CA-CFAR, 

GO-CFAR  and  SO-CFAR  detectors. 
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CHAPTER  III 


A  CELL- CENSORED  CFAR  DETECTOR  FOR 
MULTIPLE  TARGET  SITUATIONS 

3.1.  Introduction 

In  the  previous  chapter,  we  proposed  and  analyzed  the 
WCA-CFAR  detector,  which  was  shown  to  perforin  better  than  the 
CA-CFAR,  GO-CFAR  and  SO-CFAR  detectors  for  multiple  target 
situations.  In  this  chapter,  we  propose  another  scheme  for 
adaptive  cell-averaging  CFAR  detection  for  multiple  target 
situations,  where  we  censor  those  cells  which  may  contain 
interference  or  noise  spikes  in  them.  In  Section  3-2,  we 
formulate  the  problem  and  study  the  performance  of  the  pro¬ 
posed  CFAR  detector.  In  Section  3.3,  we  present  a  summary 
along  with  a  discussion. 

3.2.  Cell-Censored  CFAR  Detection 

The  cell - censored  mean-level  detector  considered  in  this 
chapter  is  shown  in  Figure  3-1.  In  this  detector,  the  idea  is 
to  delete  those  cells,  which  may  contain  interference, 
while  estimating  the  background  noise.  The  output  of  each 
range  cell  is  compared  to  a  predetermined  fixed  threshold.  A, 
to  determine  the  presence  or  absence  of  interference  in  that 


Square 

Law 

Detector 


range  cell.  The  desired  probability  of  false-alarm  at  the 
range  cells  determines  the  fixed  threshold  A.  The  comparison 
of  each  range  cell  to  A  yields  a  decision,  D^,  i=l,2,...,N, 
where  e  {0,1}.  When  D^,  i=l,2,...,N,  is  one,  it  indicates 
the  presence  of  interference  in  the  ith  cell.  This  cell  is 
then  censored  and  is  not  used  while  computing  the  estimate  of 
the  background  noise.  However,  when  D^,  i=l,2,...,N,  is  zero, 
it  indicates  the  absence  of  interference  in  the  ith  cell  and 
this  cell  is  used  in  the  estimation  of  the  background.  The 
estimate  of  the  background  noise  is  used  to  set  the  adaptive 
threshold,  Q.  As  before,  Q  is  scaled  by  a  threshold  multi¬ 
plier,  T,  prior  to  comparison  with  the  output  of  the  center 
tap  to  yield  the  desired  overall  probability  of  false-alarm 
at  the  output.  As  seen  in  the  previous  chapters,  presence  of 
interference  in  the  range  cells  raises  the  threshold  unneces¬ 
sarily  and  consequently,  it  reduces  the  probability  of 
detection  and  the  probability  of  false-alarm  drastically. 

The  proposed  cell -censoring  scheme  prevents  the  adaptive 
threshold  from  becoming  too  high  by  eliminating  the  cells 
which  may  contain  interference  and  thereby  improves  the 
system  performance. 

The  leading  range  cells,  for  which  an  absence  of 
interference  is  decided,  form  the  leading  reference  win¬ 
dow  U,  while  the  lagging  range  cells,  for  which  an  absence 
of  interference  is  decided,  form  the  lagging  reference 
window  V,  i . e . , 
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(3-1) 


I  (1  -  D  )  q 

i=  (N/  2)  + 1  1  1 


(3-2) 


Clearly,  when  is  one,  the  corresponding  q^  drops  out  of 
the  summation.  As  in  the  conventional  cell-averaging  CFAR 
detector,  the  adaptive  threshold,  Q,  is  obtained  by  summing 


the  reference  windows  U  and  V.  Q  is  then  scaled  by  the 


threshold  multiplier,  T,  in  order  to  achieve  CFAR  at  the  de¬ 


sired  value.  The  output,  qQ,  from  the  cell  under  test  (center 


tap)  is  compared  to  the  threshold,  TQ,  to  yield  the  final 


decision  D  ,  D  e  {0,1} 
o  o 


We  assume  that  the  target  to  be  detected,  or  the  primary 


target  at  the  test  cell,  is  slowly  fluctuating  in  accordance 


with  Swerling  target  model  I.  We  also  assume  that  the  inter¬ 


fering  targets  are  of  Swerling  type  I.  Then,  the  probability 


density  function  of  the  normalized  output  of  a  range  cell 


without  interference  is 


PQ.  (c*i^  =  e  1  »  -  0  » 


(3-3) 


i=l, 2 , . . . ,N 


and  the  probability  density  function  of  the  normalized  output 
of  a  range  cell  with  interference  is  given  by 
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i=l , 2 ,  .  .  .  ,N 


where  1^,  i=l , 2 , . . . ,N ,  denotes  the  SRN  of  the  interference 

(if  present)  in  the  ith  range  cell  and  the  superscript  I  in 

(°h)  indicates  the  presence  of  interference  in  the  ith 
yi  1 

range  cell. 

Next,  we  evaluate  the  performance  of  the  proposed 
system  when  there  might  be  only  one  interfering  target  with 
an  SNR  I  in  one  of  the  taps  of  the  leading  reference  window  U. 
Without  loss  of  generality,  we  assume  that  it  is  in  the  first 
cell.  Then  U  and  V  are  expressed  as 
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U  =  y  q.  +  2 


(3-5) 
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(3-6) 


where  Z  corresponds  to  the  output  of  the  first  cell  which 


may  contain  interference,  i.e., 


Z  *  (1  -  D:)  q: 


(3-7) 


The  discrete  probability  density  function  of  the  random 
variable  is  given  by 


PD1Cdl)  =  P1  6(dl  ‘  1)  +  (1  -  Pi)  St^) 


(3-8) 
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where,  6(x)  is  the  Kronecker  delta  function  such  that, 


<$(x  -  x  )  = 


1,  if  x  *  x 
0,  if  x  f  x 


(3-9) 


is  the  probability  that  the  decision  is  one,  and 
c 

(1  -  P^)  the  probability  that  is  zero.  The  superscript  c 
in  P^  indicates  that  it  is  the  probability  associated  with  an 
individual  cell.  We  may  express  P^  as 

P^  =  Pr(D1  =  1|hJ)  P(H*)  +  Pr(D1  =  1|hJ)  P(hJ) 


=  P^  P?  +  P^  P1 
1  Fi  o 


(3-10) 


where  hr  is  the  hypothesis  that  no  interference  is  present 
in  the  first  cell  while  the  hypothesis  H*  indicates  that  in¬ 
terference  is  present  in  the  first  cell.  Pp  is  the  pre¬ 
specified  probability  of  false-alarm  at  the  first  cell  which 
determines  the  fixed  threshold  A.  P1  and  P?  are  the  prior 
probabilities  indicating  the  absence  and  presence  of  the 
interfering  target  of  SNR  I  at  the  first  cell,  respectively. 
The  random  variable  Z  of  equation  (3-7)  is  defined  as 


Z  =  GQ1 


where 


G  =  1  -  Dx,  G  £  {0,1} 


(3-11) 


(3-12) 


If  is  one,  G  is  zero  and  thereby  implying  that  Z  is  zero. 
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then  from  Appendix  A, 
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Pn  ~  jj  ui  e  i  ui  1  ^ 

ui  1  r(£  -l)  1  1 


(3-18) 


The  probability  density  function  of  V  is  as  defined  in  equa¬ 
tion  (1-8).  The  adaptive  threshold  Q  is  the  sum  of  ,  V 
and  Z ,  i . e . , 


Q  =  U,  +  V  +  Z 


(3-19) 


Then,  the  probability  density  function  of  the  adaptive 
threshold,  Q,  is  the  convolution  of  the  probability  density 
functions  of  Up  V  and  Z,  i.e.. 


pq (q)  =  Py  (ux)  *  Py(v)  *  Pz(z) 


(3-20) 


where  *  denotes  convolution  as  before.  If  we  define  R  as 


R  -  U,  +  V  , 


(3-21) 


then  the  probability  density  function  of  R  is  (see  Appendix  A) 


pR(r)  =  ntr-n  rN'2  e’r  -  r  >  0 


(3-22) 
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The  probability  density  function  of  the  adaptive  threshold,  Q, 
becomes 


j  PR^r 5  PzCq-r)  dr 


(3-23) 


Substituting  equations  (3-14)  and  (3-22)  into  equation  (3-23) 
and  solving  the  integral,  we  obtain 


PQ(q) 


rl  N-2 

TTN-T)  q  e 


1-P 

■q  .  1  N-l 

rW  q 

q  >  0 


(3-24) 


To  obtain  the  expression  for  the  probability  of  detection, 
we  substitute  equation  (3-24)  into  equation  (1-4)  and  solve 
the  integral  to  obtain  the  probability  of  detection,  PD,  to 


pd  "  PJ  +  C1  ‘  pi>  f3-25^ 


The  probability  of  false-alarm  is 


PF 


(1+T) 


1  -  P. 


(1  +  T) 


(3-26) 


To  evaluate  the  performance  of  the  proposed  cell- 
censored  CFAR  detector,  we  obtain  some  numerical  results. 

We  assume  that  the  desired  overall  probability  of  false-alarm, 
Pp,  is  10  the  number  of  range  cells,  N,  is  4;  while  the 
SNR  of  the  interfering  target  is  one  tenth  of  the  SNR  of  the 
primary  target.  From  Figure  3' 2,  we  observe  that  as  \  is 
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increased  (Pp  is  reduced),  the  probability  of  detection,  PD, 
improves.  We  also  notice  from  Figures  3-3,  3-4  and  3-5  that 


the  performance  of  the  conventional  cell -averaging  CFAR  de¬ 


tector  is  better  than  the  performance  of  the  cell -censored 


CFAR  detector  for  small  values  of  the  target  SNR.  As  the  SNR 


of  the  interfering  target  increases,  at  a  certain  point  the 


cell  with  the  interfering  target  is  censored  and  the  perform¬ 


ance  of  the  cell -censored  CFAR  detector  becomes  better  than 


that  of  the  CA-CFAR  detector.  The  point  at  which  this  cross¬ 


over  takes  place  depends  upon  the  value  of  Pp  .  It  should  be 

c  1 

noted  that  if  we  make  Pp  smaller,  the  crossover  point  would 


occur  at  a  smaller  value  of  the  SNR  of  the  interfering  target 


3.3.  Summary  and  Conclusions 


In  this  chapter,  we  have  proposed  a  CFAR  detector  using 
a  cell-censoring  scheme  for  multiple  target  situations.  We 


analyzed  the  system  for  the  case  when  no  more  than  one  inter¬ 


fering  target  may  be  present.  The  analysis  for  situations 
when  more  than  one  interfering  targets  are  present  can  be 


carried  out  in  a  similar  manner.  Numerical  results  indicate 


that  the  proposed  scheme  is  effective  in  a  multiple  target 


environment. 
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Fig.  3-3.  A  Performance  Comparison  of  the  Cell-Censored  (CCJ 
CFAR  Detector  and  the  CA-CFAR  Detector. 

PF=10‘5,  Pp  =10'6,  P0  =  0.5  and  I  =  S  i  10. 
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CHAPTER  IV 


CA-CFAR  DETECTION  WITH  DISTRIBUTED  RADARS 
AND  DATA  FUSION 

4.1.  Introduction 

As  indicated  in  Chapter  One,  a  lot  of  work  on  distri¬ 
buted  detection  has  been  reported  in  the  literature  [18-30]. 
Also,  an  extensive  amount  of  research  has  been  performed  on 
cell-averaging  CFAR  detection  using  a  single  sensor  [1-17]. 
To  our  knowledge,  no  work  is  reported  on  CFAR  detection 
using  multiple  sensors  and  data  fusion.  The  goal  of  this 
chapter  is  to  develop  the  theory  of  cell -averaging  CFAR 
detection  using  multiple  sensors  and  data  fusion.  In  Sec¬ 
tion  4.2,  we  formulate  the  problem  of  distributed  CFAR 
detection  with  data  fusion.  In  Section  4.3,  we  shall  assume 
that  the  fusion  rules  at  the  fusion  center  are  known.  Spe¬ 
cifically,  we  consider  the  "AND”  and  the  "OR"  fusion  rules 
at  the  data  fusion  center.  We  obtain  the  optimum  threshold 
multipliers  of  the  individual  detectors  and  derive  an  expres 
sion  for  the  probability  of  detection  at  the  data  fusion 
center  for  the  given  fusion  rules.  As  an  illustration,  we 
study  the  performance  of  a  distributed  CFAR  detection  system 
with  two  detectors  and  data  fusion.  Numerical  results  show¬ 
ing  the  improvement  of  the  performance  for  the  distributed 
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multiple  sensor  system  with  data  fusion  over  a  single  sensor 
system  are  also  presented.  In  Section  4.4,  we  consider  the 
optimum  design  of  the  overall  system,  i.e.,  we  obtain  the 
optimum  threshold  multipliers  at  the  individual  detectors  as 
well  as  the  optimum  fusion  rule  at  the  data  fusion  center.  An 
example  is  also  presented  for  illustration.  In  Section  4.S, 


we  present  a  summary  along  with  a  discussion. 

4.2.  Distributed  CFAR  Detection  with  Data  Fusion 

We  consider  n  distributed  CA-CFAR  detectors  with  a 

data  fusion  center  as  shown  in  Figure  4-1.  It  is  assumed  that 

the  number  of  range  cells  at  the  ith  detector  is  N.  ,  i=l,2, 

...,n.  The  target  to  be  detected  is  a  slowly  fluctuating 

target  model  of  Swerling  type  I.  The  target  is  embedded 

in  a  white  Gaussian  noise  of  unknown  level.  Let  the 

probability  of  false-alarm  and  the  probability  of  detection 

at  the  individual  detectors  be  denoted  by  P„  and  , 

Fi  Di  2 

i=l,2,...,n,  respectively.  If  the  average  noise  power  is  a  , 
then  the  conditional  probability  density  function  of  the  test 
statistic  q*  from  the  test  cell  of  detector  i , i=l , 2 , . . . ,n, 
is  given  by 


A  °2 


CA-CFAR 

Cbn 


m 
m 


Data  Fusion  Center 


Fig.  4-1.  Distributed  CA-CFAR  Detection  with 
Data  Fusion. 
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where  S i,  i=l,2,...,n,  is  the  target  SNR  at  each  CA-CFAR 

detector.  The  hypothesis  Hq  represents  the  case  of  noise 

alone,  while  hypothesis  H.^  represents  the  noise  plus  target 

signal  case.  To  simplify  the  mathematical  derivations,  we 

assume  S^  =  S2  =  ...  =  Sn  =  S,  where  S  is  the  target  SNR. 

Results  for  the  case  of  unequal  target  SNR's  can  be  obtained 

in  a  straightforward  manner.  The  probability  of  detection, 

Pn  ,  for  deteitor  i,  i=l,2,...,n,  is  given  by 
Di 


Pr(Qi  >  T.q^Q1,  P  ^q1)  dq1  (4-2) 


where  is  the  scaling  factor  at  the  CA-CFAR  detector 
i,  i  =  l,2,...,n,  and  P  .  (q1)  denotes  the  probability  density 

Q1 

function  of  the  adaptive  threshold  at  the  ith  CA-CFAR  de¬ 
tector.  Also, 

00  •  •  t  q1 

PrCQ(j  >  T.qilQ1,  Hx)  =  f  P  i  (qjlHi)  dqj  =  exp  [-  ^-] 

T  q1  VH1 

(4-3) 

Since  the  noise  samples,  for  each  CA-CFAR  detector,  are 
identically  distributed,  the  probability  of  detection  of  the 
individual  detectors  can  be  written  as  (see  Appendix  A  for 


details) 


N. 

_  (1+S)  1 

.  W7 

1  (1+S+Ti)  1 


,  i=l,2 . n 


(4-4) 
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are  denoted  by 


D. 

1 


In  order  to  be 


detection,  P^, 


Then,  we  may  express  PQ,  P^,  and  Pp  as  follows 


PM  l  P0D  MD 


PF  ^  P1DFD 


and , 


P  =  1  -  P 


(4-7) 


(4-8) 


(4-9) 


where , 

7  =  summation  over  all  possible  values  of  D. 

D 

The  transition  probabilities  Pqq  and  P^  are  determined  by 
the  given  fusion  rule.  Since  D  can  take  2n  possible  values, 
there  are  2n  possibilities  for  Pqd  and  P.^.  The  goal  is  to 
maximize  the  overall  probability  of  detection  while  keeping 
the  overall  probability  of  false-alarm  constant.  To  do  this, 
we  use  the  calculus  of  extrema  and  form  the  objective  function 


’V  - 


PD(S’T1’ 


’V 


K [Pp(T1,T2, 


•V 


(4-10) 


where  v  is  the  desired  false-alarm  probability  at  the  data 
fusion  center,  £  is  the  Lagrange  multiplier,  and  T^, 
i=l,2,...,n,  are  the  threshold  multipliers  at  each  detector. 
To  maximize  Pp(S,T^ , . . . ,T  )  ,  subject  to  the  constraint  that 
Pp(Tp,T2, . . . >Tn)  is  a  constant,  we  must  maximize  the  objec¬ 
tive  function  J (Tj ,T2 , . . . ,T  ) .  We  set  the  derivative  of 
J(Tl*T2» • • • »Tr)  with  respect  to  T^,  i=l, 2, . . . ,nf equal  to 
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zero  and  solve  the  following  system  of  n  nonlinear  equations 
in  n  unknowns . 


3J(T1,T2,...,Tn) 


TTT 

J 


=  0  ,  j  =  1 , 2, . . .  ,n 


(4-11) 


Once  the  threshold  multipliers,  T\,  i=l,2,...,n,  are  obtained, 
all  the  Pp  's  are  fixed  and  the  optimum  results.  Now,  we 

give  specific  results  for  the  "AND"  and  the  "OR"  fusion  rules. 
We  also  find  the  optimum  threshold  multipliers  so  as  to  maxi¬ 
mize  Pn  while  P„  is  maintained  at  the  desired  value. 

D  r 


4.2-1  "AND"  Fusion  Rule 


In  Table  4-1,  we  present  the  "AND"  fusion  rule.  From 
this  table,  we  see  that  the  global  decision  at  the  data 
fusion  center  is  one  only  if  all  of  the  detectors  decide  a 
one.  The  transition  probabilities  are 


and 


if  D  -  [1,1, — ,1]T 
otherwise 


|  1  ,  if  D  -  [1,1,...,1]T 
0  ,  otherwise  . 


(4-12. a) 


(4-12. b) 


Substituting  equations  (4-12. a),  (4-12. b),  and  (4-9)  into 
equations  (4-7)  and  (4-8)  and  rearranging  terms,  and  Pp 
can  be  written  as 
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. 


That  is, 


=  n 


(4-1. 


(1  +  S+T^  x 


Pp  =  .  - 

F  i  =  l 


C  4  - 1 1 


(1+Ti)  1 


Substituting  equations  (4-15)  and  (4-16)  into  equation 
(4-10),  the  objective  function  is, 


J(T1,T2,...,Tn)=  ^  ^  NV 

(1+S+Ti)  1 


+  - NT 

(1+Ti)  1 


(4-1 


Taking  the  derivative  of  J  (T-^  ,^2 ,  •  •  •  ,Tn)  with  respect  to 
Tj ,  j*l,2, . . . ,n,  and  setting  it  equal  to  zero,  we  get 


Fig.  4-2.  Distributed  CA-CFAR  Detection 

with  Data  Fusion  Using  Two 
Detectors . 


PF  =  P100^1“PF1^-1'PF21  +  P101  ^•1'PF.  PF„  +  P110PF,  1PF_ 
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’AND"  Fusion  Rule 


From  the  "AND"  fusion  rule  at  the  data  fusion  center, 


the  transition  probabilities  are 


P110 


(4-29. a) 


(4-29. b) 


Thus,  the  overall  probability  of  detection,  P^,  and  the 
overall  false-alarm  probability,  Pp,  become 


PD  -  PDl  PDZ 


P  =  P  P 
F  *F  *F 


(4-30) 


(4-31) 


The  objective  function  is 


J(TlfT2)  = 


(1  +  S) 


Nl  +  N2 


IN  1  IN  «  IN  ,  IN  y 

(1  +  S  +  Tp)  i(l  +  S  +  T2)  Z  ( 1  +  T  x )  A(l+T2) 


(4-32) 


where,  and  N2  are  the  number  of  reference  cells  for  de¬ 
tector  one  and  detector  two  respectively.  We  maximize 
^^l’P2^  with  respect  to  T^  and  T?  by  setting  the  derivatives 
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equal  to  zero,  i  .  e . , 


3JCT1,T2)  3J(T1,T2) 


(4-33) 


The  resulting  equations  are 
N  +N 

(1  +  S)  ^  r  1  „  fA  ... 

R7+T  N“  S  nT^I  R7  '  0  (4'34) 

(1+S+T1)  1  (1+S+T? )  (1+T  )  1  (1+T2) 


(1  +  S) 


Nl  +  N2 


fj  “  M  +f  ^  M  KT 

(1+S+Tj)  1(1+S+T2)  2  (1+T1)  1(1+T2)  2 


living  for  and  T2  subject  to  the  constraint 


1  1  "? 
(1  +  T1)  d+T2)  Z 


7X  =  0  (4-35) 


(4-36) 


we  obtain, 


T1  =  T2  =  -!  +  v 


Nl  +  N2 


(4-37) 


The  performance  of  the  system  in  terms  of  Pn  versus  the 


target  SNR  for  =  4,  N2  =  6,  and  =  10  is  plotted  in 
Figure  4-3. 


'OR”  Fusion  Rule 


In  this  case,  the  transition  probabilities  are 


Using  the  constraint  that 


pf  =  — ~ — ?r  * - ^tt  - - r - tt*  v  t4'43’ 

(l+Tj)  1  fl*T2)  2  Cl+Tj)  1(l+T2)  2 
we  obtain  in  terms  of  as 

N,  1 
(1+T  )  z  -  1  W7 

T,  =  -  1  +  [ - - ]  1  (4-44) 

v  (1+T2)  2  -  1 

Substituting  equation  (4-44)  into  equation  (4-43)  and  using 
an  iterative  procedure,  we  obtain  a  value  for  and  then 
For  =  4,  =  6,  and  Pp  =  10  ^ ,  the  values  of  and 

are  11.597  and  4.050,  respectively.  The  numerical  results 
showing  the  performance  are  plotted  in  Figure  4-4.  In  the 
next  subsection,  we  comment  on  the  performance  of  the  "AND" 
and  the  "OR"  fusion  rules. 

4.2-4.  A  Comment  on  the  Performance  of  the 
"AND  and  "OR"  Fusion  Rules 

In  Figure  4  ^ ;  we  plot  the  system  performance  for  the 
"AND"  and  the  "OR"  fusion  rules.  As  evident  from  the  figure 
when  the  target  SNR  exceeds  14.5  dB  (approximately),  the 
performance  of  the  "OR"  fusion  rule  becomes  better  than  the 
performance  of  the  "AND"  fusion  rule.  This  is  easily  demon¬ 
strated  analytically.  Let  the  overall  probabilities  of 
detection  at  the  data  fusion  center  for  the  "AND"  and  the 
"OR"  fusion  rules  be  denoted  by  Pn  and  Pn,  respectively. 


>y 


*D  -  D 


(4-45) 


if  and  only  if, 
Nn 

(1+S) 


N. 


"IT 


(l+S) 


it: 


(l  +  S) 


Nl  +  N2 


TT 


(1+S+Tj) 


(1+S+T2) 


(1+S+T1)  1(1+S+T2)  2 


73T  - 


(l+S) 

- n 

(l+S+T1) 


Nl  +  N2 


1(i+s+t2) 


TT 


(4-46) 


Solving  this  nonlinear  equation,  we  obtain  S=14.57  which 
agrees  with  the  graphical  result. 


4.3.  Optimization  of  the  Overall  System 

In  the  previous  section,  we  studied  the  problem  of 
CA-CFAR  detection  with  multiple  sensors  and  data  fusion. 

Given  the  fusion  rule  at  the  data  fusion  center,  we  maximized 
the  overall  probability  of  detection  by  obtaining  the  opti¬ 
mum  threshold  multipliers  at  the  individual  CA-CFAR  detectors, 
while  the  overall  probability  of  false-alarm  was  maintained 
at  the  desired  value.  In  this  section,  we  will  also  consider 
the  design  of  the  optimum  fusion  rule.  Again,  the  goal  will 
be  to  maximize  the  overall  probability  of  detection  while  the 
overall  probability  of  false-alarm  is  kept  at  a  desired  value. 
While  maximizing  the  overall  probability  of  detection,  we 
will  derive  the  optimum  fusion  rule  at  the  data  fusion  center 
as  well  as  the  optimum  threshold  multipliers  at  the  indivi¬ 
dual  CA-CFAR  detectors.  We  will  use  an  approach  similar  to 
the  one  in  [35]. 

The  decision  output,  Dq ,  at  the  data  fusion  center  is  a 
function  of  the  incoming  decisions  ,  i=l,2,...,n,  i.e., 
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Dq  =  f (Dj ,D2 ,  .  .  . ,Dn) 


(4-47. a) 


Dq  e  {0,  1}  (4 -47 . b) 

For  ease  in  handling,  we  use  the  following  notation  for 
quantities  defined  in  equations  (4-6. b),  (4-6. c)  and 
(4-6. d), 


P1D  =  Pr<Do  = 


PD1  =  MD 


(4-48. a) 
(4-48. b) 


P  =  F 
*DO  D 


(4-48. c) 


Then,  the  overall  probability  of  detection,  Pp,  and  the 
overall  probability  of  false-alarm,  Pp,  at  the  output  of  the 
data  fusion  center  can  be  written  as 


=  l  P1D  PD1 

(4-49) 

=  d  ?1-  P-° 

(4-50) 

The  optimization  process  for  the  overall  system  will 
be  carried  out  in  two  steps.  First,  we  optimize  the  fusion 
rule  at  the  data  fusion  center  assuming  that  the  threshold 
multipliers  are  known.  This  yields  a  set  of  equations.  Then, 
we  obtain  the  optimum  threshold  multipliers  assuming  that 
the  optimum  fusion  rule  is  known,  which  gives  another  set  of 
equations  to  be  solved.  These  two  sets  of  nonlinear  equa¬ 
tions  are  then  solved  simultaneously  to  yield  the  optimum 
threshold  multipliers  and  the  optimum  fusion  rule. 
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Our  goal  is  to  maximize  J.  This  maximization  is  equivalent 
to  maximizing  P1D*  PDn  subject  to  the  constraint  P1D*PD*0S 
Dividing  equation  (4-54)  by  PD*i  we  obtain 


where 


J*  =  Pln*  +  S*  {Pln*  -  vv*}  +  KK(D*) 


(4-55) 


KK(D*)  = 


K(D*) 


Since  KK(D*)  does  not  depend  on  P1D*,  equation  (4-55)  repre 

sents  the  equation  of  a  straight  line  in  terms  of  P1T^*.  Thi 
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as  in  [35]  ,  we  obtain  the  conditions  which  the  transition 
probabilities  must  satisfy  to  yield  the  optimum  fusion  rule 


=  min  {  vv*,  1 } 


(4-56) 


ID*  =  0 


The  next  step  is  to  find  the  optimum  threshold  multi¬ 
pliers,  T^,  i=l,2,...,n.  The  probabilities  of  detection, 

PD  ,  and  the  probabilities  of  false-alarm,  Pp  ,  i=l,2,...,n, 
at  each  CA-CFAR  detector  are  functions  of  the  threshold 


multipliers,  T\  ,  i=l,2,...,n.  Therefore,  maximizing  P^, 
while  Pp  is  fixed  at  the  desired  value  v,  is  equivalent 
to  maximizing  J  with  respect  to  the  threshold  multipliers 
T^'s.  We  take  the  derivative  of  J  with  respect  to  , 
i=l , 2 , . . . ,n,  and  set  them  equal  to  zero  to  yield  the  follow¬ 
ing  set  of  equations  which  the  optimum  T^'s  must  satisfy. 


0,  x  1 , 2  ,  .  .  .  ,n 


(4-57) 


It  should  be  emphasized  that  the  set  of  equations  (4-56) 
and  (4-57)  are  to  be  solved  simultaneously  in  order  to  ob¬ 
tain  the  optimum  threshold  multipliers  as  well  as  the  best 
fusion  rule.  The  computations  involved  are  overwhelming  and 
it  will  require  an  extensive  computer  iterative  procedure. 
This  point  will  become  evident  in  the  simple  example  we 


consider  next 


Example 


We  consider  the  system  of  two  CA- CFAR  detectors  with 
data  fusion.  For  the  sake  of  simplicity,  let  the  transition 
probabilities  of  equation  (4-26)  be  denoted  by  cu  ,  i=l,2,3,4. 


i  .  e .  , 


al  ~  P100 


a2  P101 


a3  "  P1 1 0 


a4  =  P1 1 1 


(4-58. a) 
(4-58. b) 


(4-58. c) 


(4-58. d) 


Using  the  same  notation  as  in  the  example  presented  in 
Section  4.2,  we  can  write  the  overall  probability  of  detec- 
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tion,  Pp,  and  the  overall  probability  of  false-alarm,  Pp, 


PD  “  al^1'PD1^1'PD2^  +  a2^1‘PD1^)PD2  +  ot3PD1  ^1_PD2^  +  a4PDPD2 

(4-59) 


P„  =  an  (l"Pc  )  (1  -Pr;  )  +  Ct^Cl-Pp  )P-  +  a,Pp  Cl -Pt-  )+  ct.Pn  Pc 

F  lv  Fi  F2  2  Fj-'  F 2  3  F1'-  F2  4F1F2 

(4-60) 

Assuming  that  the  target  SNR  is  the  same  at  each  CA-CFAR 
detector,  denoted  by  S,  the  probabilities  of  detection  and 
the  probabilities  of  false-alarm  at  each  CA-CFAR  detector 
are  (see  Appendix  A) , 


n  r  1  +  S  ,  1 

FD1  =  L1  +  S  +  T1J 


„  r  1  +  s  ,2 

FD2  L1  +  S  +  T2j 


1  (1  +  T,)  1 


2  w2 
(1  +  T,) 


(4-61) 


(4-62) 


(4-63) 


(4-64) 


As  discussed  earlier,  to  optimize  the  entire  system  we  need 
to  maximize  the  objective  function  defined  in  equation  (4-51) 
Substituting  equations  (4-59)  and  (4-60)  into  equation  (4-52) 
we  obtain 

J  =  ai^1'PD1'PD2+PD1PD2^  +  a2^1'PD1^)PD2  +  a  3PD ^1-PD 2^+a4PD1PD  2 

+  £{a, (1-PC  -Pc  +PC  Pc  )  +  (^-(l-Pr.  )PC  +a,Pc  (1-PE  )+0l,Pc  ' 

1  F1  F  2  F1  F2  ^  F1  F2  3  F1  F2  4  F1  F2 

(4-65) 


k  >  Vk 


Dividing  equation  (4-65)  by  (1-Pn  _Pn  +Pn  Pn  )»  we  get  a  new 

1  U2  U1  U2 


objective  function  Jp  i.e.. 


where 


=  al  +  ^1  ^al  ‘  vi  ^  +  Ki 


AD1  *  1  -  pd1'pd2  *  pd1pd2 


AF,  =  1  -  P-  -Pc  +  Pc  Pc 
1  F1  F2  F1  F2 


J  = 

J1  AD 


1 

AF, 


?1  ^  AD, 


(4-66) 


(a?(l-Pp,)P  +  a,Pp  (1-PF  )  }  +  a.Pp  Pp 

/  h.  t2  ~>  ^  rx  t2 


AF 


1 


K,  = 


a2(1'PD1^PD2  +  a3PD1(1‘PD2)  +  0t4PD1PD2 


AD, 


Notice  that  equation  (4-66)  is  the  equation  of  a  straight 
line  in  .  In  this  case,  the  extrema  are  at  the  end  points 
depending  on  whether  the  slope  is  positive  or  negative.  Thu 
at  the  data  fusion  center,  the  value  of  is  given  by 


(4-67) 


In  an  analogous  manner,  we  obtain 


«  t«,  «,»  ,>  4,i  iA't:*' »Vi:«,yjl»M» >*  w 


where , 


J11  =  PDX  +  511  {PF1  '  Vll}  +  K] 


„  Val  +  PF2(al'a2'a3  +  a4> 
?11  "  ?  a3'al  +  PD2Cal’a2'a3  +  a4^ 

al"V  +  PF2(a2'al) 

VU  '  a3-aL  +  PD2  (a1-a2-a3  +  a4) 

al  +  PD2(a2'al} 

KU  "  a 3-ai  +  PD2Ca1-a2-a3  +  a4) 


(4-72) 


In  a  similar  manner,  we  can  define  J22>  ?22>  v22  anc*  K22  Fo 
the  second  threshold  multiplier.  Taking  the  derivatives  of 
J.^  and  J22  with  respect  to  and  T2>  respectively,  and 
rearranging  terms,  we  obtain 


1  _  .  F 
3PC  ni 
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2  _ 

3Pl  ‘  ^22 
F2 


(4-73) 


(4-74) 


Substituting  equations  (4-61)  to  (4-64)  into  equations 


(4-73)  and  (4-74),  we  get 


N,  N.+l 

(1+S)  1  (1+Tj)  1 

(1  +  S  +  Tx) 

N~  N_  +  l 

(1+S)  “  (1+T.)  Z 


(1  +  S  +  T-) 


(4-75) 


(4-76) 


Solving  equations  (4-59)  to  (4-64),  (4-67)  to  (4-70), 
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(4-73)  and  (4-74)  by  an  extensive  iterative  procedure,  for 
Ni  =  6,  Nz  =  4,  ve  obtain  the  optimum  Pp  versus  Pp  curves 
shown  in  Figures  4-6  and  4-7,  for  S  =  5  dB  and  S  =  10  dB, 
respectively.  To  obtain  the  overall  probability  of  detec¬ 
tion,  PD,  for  a  specific  probability  of  false-alarm  at  the 

data  fusion  center,  we  read  it  directly  from  the  graphs. 

-  4 

From  Figure  4-6,  we  read  that  for  Pp  =  10  ,  Pp  is  .451, 

which  is  equal  to  the  Pp  obtained  from  the  "AND"  fusion 
rule . 

4 . 4  Summary  and  Discussion 

In  this  chapter,  we  have  developed  the  concept  of  dis¬ 
tributed  CA-CFAR  detection  with  data  fusion,  where  detec¬ 
tion  decisions  are  transmitted  from  each  CA-CFAR  detector 
to  the  data  fusion  center.  The  overall  decision  is  obtained 
at  the  data  fusion  center  based  on  some  "k  out  of  n"  fusion 
rule.  In  the  first  part  of  the  chapter,  we  obtaind  the 
optimum  threshold  multipliers  at  each  CA-CFAR  detector, 
given  the  fusion  rule  at  the  data  fusion  center.  In  the 
latter  part  of  the  chapter,  we  optimized  the  entire  system. 
That  is,  we  obtained  the  optimum  threshold  multipliers  at 
each  CA-CFAR  detector,  as  well  as  the  optimum  fusion  rule  at 
the  data  fusion  center.  Although  the  computations  involved 
required  an  extensive  iterative  procedure,  numerical  re¬ 
sults  were  obtained  to  show  the  performance  of  the  systems 
proposed . 


FALSE -ALARM  PF 


PD  Versus  Pp  for  the  Overall  Optimum  System 
S  =  10  dB. 


CHAPTER  V 


ADAPTIVE  CFAR  DETECTION  WITH  MULTIPLE 
BACKGROUND  ESTIMATORS 


5 . 1  Introduction 

In  the  last  chapter,  we  studied  the  use  of  multiple 
sensors  with  data  fusion  for  adaptive  CFAR  detection.  De¬ 
cisions  of  the  individual  CFAR  detectors  were  transmitted  to 
a  data  fusion  center  where  they  were  combined  according  to 
some  fusion  rule  to  yield  a  final  decision.  In  this  chapter, 
we  will  still  use  multiple  sensors  as  shown  in  Figure  5-1. 
However,  rather  than  the  detector  decisions,  the  individual 
sensors  compute  and  transmit  estimates  of  the  level  of  the 
background  noise  to  the  CFAR  detector.  These  estimates  are 
used  to  set  the  adaptive  threshold  required  for  the  desired 
probability  of  false-alarm  at  the  CFAR  detector.  In  Section 
5.2,  we  formulate  the  problem,  discuss  the  operation  of  the 
system  and  derive  expressions  of  the  probability  of  detection 
for  CA-CFAR,  GO-CFAR,  and  SO-CFAR  detection.  An  example  is 
presented  in  Section  5.3.  In  Section  5.4,  we  suggest  a 
hybrid  system  which  combines  the  concept  developed  in  Chapter 
four  and  the  concept  used  here  in  Section  5.2. 


Phenomenon  H 


Fig.  5-1.  CFAR  Processor  with  Multiple  Background 
Estimators . 
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5.2  Adaptive  CFAR  Detection  with  Multiple  Background 
Estimators 

A  more  detailed  block  diagram  of  the  system  under 
consideration  is  shown  in  Figure  5_2.  It  consists  of  (n+1) 
sensors,  n  of  which  are  used  solely  for  the  estimation  of 
the  unknown  level  of  the  white  Gaussian  background  noise. 
These  sensors  are  denoted  by  E^ ,^2 > • • • ,En-  The  sensor  E^, 
i=l,2,...,n,  obtains  the  maximum  likelihood  estimate  (MLE) t 
,  i=l,2,...,n,  of  the  unknown  level  of  noise  (we  refer  to 
VL  as  the  MLE  even  though  strictly  speaking  it  is  times 
the  maximum  likelihood  estimate,  where  is  the  number  of 
range  cells  at  the  ith  estimator).  As  indicated  in  Figure 
5-2,  the  center  tap  and  the  guard  cells  are  not  used  during 
the  computation  of  at  the  sensor  E^.  After  this  data 
compression,  the  MLE  ,  i=l,2,...,n,  is  transmitted  to  the 
CFAR  detector  CD.  CD  computes  its  own  ML  estimate  of  the 
background  noise  and  obtains  the  reference  windows  U  and  V. 
The  MLE's  W. ,W- , . . . ,W  are  combined  with  the  reference  win- 
dows  U  and  V  at  sensor  CD  where  a  decision  about  the  target 
is  to  be  made.  The  MLE's  U,V  and  ,  i=l ,  2  , . . . ,n,  are  com¬ 
bined  according  to  the  desired  selection  logic  CA,  GO  or  SO 
to  yield  the  adaptive  threshold  Q.  Q  is  then  scaled  by  the 
threshold  multiplier,  T^,  so  that  the  probability  of  false- 
alarm  at  the  output  of  the  CFAR  detector  at  CD  is  maintained 
at  the  desired  value.  The  test  statistic,  q  ,  from  the  cell 
under  test  at  the  CFAR  detector  is  compared  to  the  threshold 
T^Q  in  order  to  make  a  decision. 


The  target  to  be  detected  is  assumed  to  be  of  Swerling 
type  I  with  an  SNR  S  embedded  in  a  white  Gaussian  noise.  The 
noise  samples  of  the  reference  cells  of  the  sensors  are 
assumed  to  be  independent,  identically  distributed  with  the 
probability  density  function  defined  as  in  equation  (1-5). 

The  probability  density  function  and  the  cumulative  distribu¬ 
tion  function  of  the  estimates  W^,  i=l  ,2  , .  .  .  ,n,  are  (see 
Appendix  A  for  details)  . 
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(5-1) 


(wA) 
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j  =  0 


L  0, 

i*— 1,2,.  . .  ,n 


(5-2) 


where  N^ ,=1 , 2 , . . . ,n ,  denotes  the  number  of  reference  cells 

at  the  sensor  E^ ,  i=l,2,...,n.  The  reference  windows  U  and  V 

at  CD  are  assumed  to  have  the  same  number  of  range  cells 
N 

equal  to  -j.  For  ease  of  notation  in  deriving  the  adaptive 
threshold,  we  let 


*  =  N  (5-3) 

2  o  v  1 

Now,  we  obtain  the  probability  of  detection  for  the 

CFAR  processor  with  multiple  background  estimators  shown  in 

Figure  5-2.  The  cell -averaging  selection  logics  considered 

are  the  CA,  GO  and  SO. 


5.2-1  CA-CFAR  Detection 


In  this  case,  the  statistic  Q  is  the  sum  of  the  MLE's 


(5-4) 


Q  =  U  +  V  +  W1  +W  + 


where  U,  V,  and  W^,  i  =  l,2,...,n,  are  statistically  inde¬ 
pendent  random  variables.  Let  us  define 


W  =  U  +  V 
o 


(5-5) 


Then,  the  probability  density  function  of  WQ  is  the  convolu¬ 
tion  of  the  probability  density  functions  of  U  and  V.  Using 
equations  (1-7)  and  (1-8) ,  and  rewriting  the  convolution  to 
obtain  the  probability  density  function  of  W  ,  we  have 

-w  wo  -  -  1  N  _  i 

pw  C«o>  '  rw rrtUT  e  °  I  “2  (“'"s’2  du 

O  v  O  O 


(5-6) 


Solving  the  integral  and  rearranging  terms,  we  obtain 
PW0fwo:)  -  TTi^T  e  °  "o°  ■  "o  ^  0 


(5-7) 


The  probability  density  function  of  Q,  Pq(q),  is  given  by 


the  following  n-fold  convolution. 


PQ  Cq)  '  PW  ("o’  *  PW, ("l’  *  • ' '  *  PW  ("„’  1 

o  1  n 

For  the  sake  of  clarity,  we  define  the  following  random 
variables . 


Qi  =  wo  +  wi 


(5-8) 


Qi  =  Qj_i  +  »  i=2 ,3 , . . . ,n-l 


Q  =  Q  .  +  W 
x  xn-l  n 


(5-9) 
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The  probability  density  function  of  Qj  is  given  by  the 


following  convolution 


PQ  ^ 


PW  tw0>  PW,tqfV  dW0 

o  1 


(5-10) 


Substituting  equations  (5-1)  and  (5-7)  into  equation  (5-10), 


we  obtain 


PQ1Cqi)  “  r(Nx)  r(No; 


Solving  the  integral,  we  get 


■N  -1 


(qr"0)  1 


N,  -1 


(5-11) 


N  +N..-1 

qx  0  B(No,N1) 


(5-12) 


where  B(x,y)  is  the  incomplete  Beta  function.  B(x,y)  can  be 


defined  in  terms  of  the  Gamma  function,  i.e.. 


B(x,y)  = 


IWn 

r (x+y' 


(5-13) 


Substituting  equation  (5-13)  into  equation  (5-12),  the 
probability  density  function  of  becomes 


PQ  ^ql^ 


q  i  W1 


1  +  N, 
o  1' 


»  i0 


(5-14) 


Similarly,  we  can  derive  the  probability  density  function  of 


Qt »  Pn  ( q 9 )  •  Since 


Q?  =  Q,  +  W. 


(5-15) 


and  and  W ^  are  statistically  independent. 
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iolving  the  integral  and  using  equation  (5-13),  Pn  (q_)  be- 

X?  4 


comes 


PQ2^q2^  = 


r(No+N1+N2) 


- q  _  N  +N,+N,-1 
^  2  o  1  2 

e  q  2 


(5-17) 


Continuing  in  this  manner,  the  probability  density  function 


LiA 

of  the  adaptive  threshold,  Pq  (q) ,  can  be  seen  to  be 


where , 


r,  CA  r  1  „-q  _  4+  1 

pq  (q;)  -  tt^t  q 


*  -  I  N 

i  =  0  1 


(5-18) 


Substituting  equations  (5-18)  and  (1-3)  into  equation  (1-2), 

CA 

the  probability  of  detection,  Pn  ,  is 


TCAq 

1  +  S  1  -q  ip-1  , 

•  nw e  q  dq 


(5-19) 


This  results  in 


1  +  S  ■,  ip 


(5-20) 


.CA 


The  probability  of  false-alarm,  Pp  ,  is  obtained  by  just 
setting  the  target  SNR  S  to  be  zero,  i.e., 


P_  =  - = - r 

*  tca]* 


5.2-2  GO-CFAR  Detection 


In  this  case,  the  adaptive  threshold  Q  is 
Q  =  MAX(U,V,W1,W2, . . . ,Wn) 


(5-21) 
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Since  the  random  variables  U,V,  and  W^,  i=l,2,...,n,  are 
independent  of  each  other,  the  cumulative  distribution  func¬ 
tion  of  Q  is 

Fq0 (q)  =  Pr  (Max(U,V,W1,W2,...,Wn)  <  q} 

=  Pr  (U  £  q ,  V  £  q ,  £  q ,  Vi ^  <_  q,  .  .  .  ,Vi^  <_  q} 

=  Pr  {U£q  }  *Pr  {V£q}  •Pr{W1£q  }  •  Pr  {W2_<q  } .  .  .  Pr(Wn£q} 

=  Pjj (q)  Fv(q)  Fw  Cq)  Fw  (q)  .  .  .  Fw  (q)  (5-22) 

12  n 

The  probability  density  function  and  the  cumulative  distri¬ 
bution  function  of  U  and  V  are  as  defined  in  equations  (1-7) 
to  (1-10).  We  define  the  notation 


u  = 

w 

0 

and 

7=  N 

2  o 

(5- 

•23) 

V  = 

Wn+1 

and 

N  M 

1  ~  Nn+1 

(5- 

■24) 

Using  equations  (5-23)  and  (5-24)  in  equation  (5-22),  the 
cumulative  distribution  function  becomes 

F0°(q)  =  Fw  (q)  Fw  (q)  Fw  (q)  ...  Fw  (q)  Fw  (q)  (5-25) 
v  o  1  2  n  n+1 

Then,  the  probability  density  function  of  the  adaptive  thresh 
old  is  obtained  by  taking  the  derivative  of  equation  (5-25) 


(q)  = 

pw  (q) 

0 

Vq> 

Fw  (q) 
w2 

• • •  F^  (q ) 
n 

Fw 

wn+l 

+ 

Fw  (q) 
0 

PWi(q) 

\W 

•  •  •  F^  ( q ) 
n 

Fw  (q) 

wn+l 

+ 

+ 

f^  (q ) 
0 

• 

• 

FWX  ^ 

V"1 

*  •  •  F^  (q) 
n 

Fw  (q) 
wn+l 

+ 

fw  (q) 

o 

Vq) 

fw2  ^ 

* • •  (q ) 

n 

FW  Cq) 

wn+ 1 

+ 

Fw  (q) 
0 

FWX  ^ 

F^tq) 

...  Fw  (q) 
n 

PW  (q) 

wn+l 

Now,  we  substitute  equations  (1-7)  to  (1-10),  (5-1),  and 
(5-2)  into  equation  (5-26)  to  obtain 
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are  defined  in  equations  (1-7)  and  (1-8),  and  the  probability 
density  function  of  W  is 


pw(w)  = 


N+l  -w 
w  e 


w  >  0 


(5-33) 


Using  these  probability  density  functions,  we  derive  expres¬ 
sions  of  the  probability  of  detection  for  the  three  schemes, 
"A,  GO  and  SO,  discussed  above. 


CA-CFAR 


The  adaptive  threshold  Q  is  the  sum  of  U,V,  and  W,  i.e., 


Thus , 


Q  =  U  +  V  +  W 


pq  Cq)  =  pu(u)  *  Pv(v)  *  Pw(w) 


(5-34) 


(5-35) 


CA 

Using  the  results  derived  in  the  previous  section,  (q)  is 


p^A(q)  - 


2N+1  -q  .  n 
q  e  q  >_  0 


(5-36) 


Substituting  equations  (5-36)  and  (1-3)  into  equation  (1-2), 


we  obtain  the  probability  of  detection  Pn  to  be 


r,CA  r  1  +  S  •,  2N+  2 

PD  =  +  3  V  TcaJ 


while  the  probability  of  false-alarm  is 


^  +  W 


(5-37) 


(5-38) 
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For  a  probability  of  false-alarm  of  10  and  N=4,  we  plot 
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Pp  versus  the  target  SNR  to  study  the  improvement  in  the 
performance.  In  Figure  5-4,  we  plot  the  probability  of 
detection  of  the  two-sensor  system  of  Figure  5_3,  and  the 
probability  of  detection  of  a  single  sensor  system.  The 
improvement  in  the  performance  is  clearly  noticeable. 

GO-CFAR 

In  this  case,  the  adaptive  threshold  is  eqi.al  to  the 
maximum  of  U,V,  and  W,  i.e., 

Q  =  MAX ( U , V , W )  (5-39) 

The  cumulative  distribution  of  Q  is 

FQ°(q)  =  FyCq)  Fy(q)  Fw(q)  (5-40) 

while  the  probability  density  is 


Pq  Cq)  =  Pu(q)Fv(q)Fw(q)  +  Fy(q) Py (q) Fw (q)  +  Fy (q) Fy(q) Pw (q) 

(5-41) 

Substituting  equations  (1-7)  to  (1-10),  and  (5-33)  into 
equation  (5-41),  the  probability  density  function  of  the 
adaptive  threshold  becomes 


PG0 
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Performance  of  the  CA-GFAR  Detectors  for  Pp  =  10 
and  N=4. 

DS  =  CFAR  Processor  with  a  Background  Estimator. 
SS  =  Single  Sensor  CFAR  Detector. 
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(5-42 


The  probability  of  detection  is  obtained  from  equation  (1-2) 

by  using  the  results  of  equations  (5-42)  and  (1-3).  Solving 

CO 

the  integral  and  rearranging  terms,  ,  can  be  written  as 


Perform 
and  N=4 
DS  =  CF. 
SS  5  si 


SO-CFAR 


The  threshold  in  this  case,  is  the  minimum  of  U,V, 


and  W,  i . e . , 


Q  =  MIN (U ,V,W) 


(5-44) 


The  cumulative  distribution  of  Q  is  given  by 


Fg^(q)  -  Fy(q)  +  Fy(q)  +  F^(q)  -  Fyy(q,q)  -  Fy^(q,q)  -  Fy^(q,q) 


+  FUVW ^ ^ 


The  probability  density  function  is,  therefore, 


(5-45) 


Fq0 (q )  =  2P[J(q)  [1  -  Fy  (q)  -  Fw(q)  ]  +  2Pu(q)  Fy(q)  Fw(q) 


+  pw(q)t1  -  2Fy(q)]  +  Pw(q)  Fjj  Cq )  Fy(q) 


(5-46) 


Substituting  equations  (5-46),  (1-3),  (1-7),  (1-8),  (1-9), 
(1-10)  and  (5-33)  into  equation  (1-2),  and  rearranging  terms, 
the  probability  of  detection  is  obtained  to  be 

2  W2)"1  ^  r<7  ♦  j  *  k>  r  X  +  S 

pd  ■  73T  x-o  Jo  — Tttt - 'Wl^1 


CN/ 2) - 1  (N/2) - 1  (N,f,k,1)( 

i  i  — rnn - 

j  =  0  k=0  JK- 


1  +  S  -.N+i+k  +  2 
' 


(5-47) 


In  Figure  5-6,  we  compare  the  performance  of  all  of  the 
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three  schemes  for  the  two-sensor  system  considered  in  this 
example . 

5.4.  An  Extension 

In  this  section,  the  concept  of  using  multiple  back¬ 
ground  estimators  is  employed  in  conjunction  with  the  con¬ 
cept  of  distributed  detection  with  data  fusion  developed  in 
Chapter  four.  The  hybrid  system  that  we  consider  is  shown 
in  Figure  5-7.  The  system  consists  of  n  background  estimators 
El*^2,.',,^n  anc*  a  ^FAR  processor.  Each  background  estimator 
E^,  i=l,2,...,n,  transmits  the  output  from  its  test  cell 

(center  tap)  as  well  as  the  MLE  of  the  level  of  the  background 
noise,  W^,  to  the  CFAR  processor.  The  MLE's  are  combined  to 
obtain  the  adaptive  threshold  which  is  scaled  by  the  threshold 
multiplier,  T^,  to  achieve  the  desired  value  of  the  probability 
of  the  false-alarm  at  the  system  output.  At  the  CFAR  processor 
each  test  cell  qQ^,  i=l,2,...,n,  is  compared  to  the  above 
adaptive  threshold  to  yield  a  decision  D^,  i=l,2,...,n,  where 
e  (0,  1}.  In  a  manner  analogous  to  Chapter  four,  the  de¬ 
cisions  are  combined  according  to  a  "k  out  of  n"  fusion  rule 
to  yield  the  final  decision  D  .  The  performance  of  the 
system  proposed  in  this  section  is  expected  to  be  better  than 
the  system  discussed  in  Chapter  four.  This  is  due  to  the 
additional  information  about  the  background  noise  available 
at  the  CFAR  processor  while  computing  the  adaptive  threshold 
Q.  To  illustrate  the  above  concepts,  an  example  is  presented. 


Example 


In  this  example,  we  consider  the  system  consisting  of 

two  background  estimators  and  a  CFAR  processor  as  shown  in 

Figure  5-8.  At  the  CFAR  processor,  we  use  the  CA-CFAR  scheme 

and  combine  the  individual  decisions  according  to  the  "AND” 

fusion  rule.  The  result  is  compared  to  the  corresponding 

system  discussed  in  Chapter  four,  i.e.,  CA-CFAR  detection 

with  the  "AND"  fusion  rule.  Let  and  pi}  denote  the 

ul  u2 

probabilities  of  detection  corresponding  to  the  decisions 

Dj  and  respectively,  in  the  hybrid  system.  Let  Pp  and 
^  1 
P~  be  the  probabilities  of  detection  corresponding  to  the 
u2 

decisions  D^  and  D^  for  the  system  shown  in  Figure  4-2.  Let 
}i  d 

Pp  and  P°  represent  the  overall  probabilities  of  detection 
for  the  hybrid  system  and  the  system  of  Figure  4-2,  respec¬ 
tively.  Then,  as  shown  in  Chapter  four,  the  detection  proba¬ 
bilities  can  be  written  as 


1  +  SV  Nn  +N? 

r - rH— ]  1  2 


1  +  s“  +  T 


1  +  S‘  n,+n. 
[ - TT^ - ]  1  2 


1  +  S“  +  T 


1  +  N. 

t - ^  -  j  1 


1  +  S~  +  T1 


1  +  S7  N_ 

[ - X-2 - ]  2 


1  +  s2  +  T2 


C5-48) 


(5-49) 


(5-50) 


(5-51) 
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where  and  are  the  target  SNR's  at  the  two  sensors  of 
the  hybrid  system;  Sd  and  Sd  are  the  target  SNR's  at  the 
system  of  Figure  4-2;  and  are  the  threshold  multi¬ 
pliers  of  the  system  shown  in  Figure  4-2;  and  T  is  the 
threshold  multiplier  of  the  hybrid  system.  Substituting 
equations  (5-48),  (5-49),  (5-50)  and  (5-51)  into  equation 
(4-13),  the  overall  probabilities  of  detection  P^  and  Pd  can 
be  written  as 


,  1  +  S1?  1  -►  S9 

Ph  =  {  1  .  2 

^  1  +  S^  +  T  i  +  +  T 


ni+n2 


(5-52) 


j  1  +  s, 

pd  =  r _  1 

D  1  +  S  j  +  T  j 


1  +  S,  N, 

_ t - 1  L 

3  J 


1  +  S2  +  T2 


(5-53) 


In  order  to  compare  the  two  systems  under  the  same  conditions, 
we  let  S^1  =  -  Sd  «  Sd  «  S.  Then,  Pp  and  Pd  become 


Ph  =  r 1  *  -s i 

D  ll  +  S  +  TJ 


2(Nx+N2) 


Pd  = 
D 


(1  +  S) 


n1+n2 


1  1  2 
(1  +  S  +  T1)  1  (1  +  S+T2)  Z 


(5-54) 


(5-55) 


The  probabilities  of  false-alarm  are 


rh  _  1 

r  2'(N7+N  ' 

(1  +  T)  1  L 


pd  _  _ 1  _ 
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Using  Nj  -  6,  N2  =  4,  and  Pp  =  Pp  -  Pp  -  10‘4,  we 
plot  P^  and  Pp  versus  the  target  SNR  S  in  Figure  5-9.  As 
expected,  the  hybrid  system  performs  better. 

5.5.  Summary  and  Discussion 

In  this  chapter,  we  have  developed  the  concept  of  using 
multiple  background  estimators  for  CFAR  processing.  Since 
the  multiple  background  estimators  provide  more  information 
about  the  level  of  noise  in  the  background,  the  performance 
improved  as  expected.  The  model  considered  in  this 
chapter  is  applicable  to  systems  using  frequency 
diversity.  The  model  is  also  applicable  to  space 
diversity  systems  if  we  make  idealized  assumptions, 
i.e.,  perfect  synchronization  etc.  We  also  discussed 


an  extension  of  the  work  in  Chapter  four  using  multiple 
background  extimators. 


rformance 

timators 


CHAPTER  VI 


ADAPTIVE  CFAR  DETECTION  WITH  DIFFERENT 
NETWORK  TOPOLOGIES 

6.1.  Introduction 

In  the  last  two  chapters,  we  have  considered  the  problem 
of  multisensor  CFAR  detection  with  data  fusion.  In  Chapter 
four,  the  decisions  from  distributed  CFAR  processors  were  com¬ 
bined  at  the  data  fusion  center  according  to  a  "k  out  of  n" 
fusion  rule  to  yield  the  overall  decision.  In  Chapter  five, 
multiple  sensors  transmitted  the  estimates  of  the  background 
noise  to  the  CFAR  processor  where  an  adaptive  threshold  based 
on  these  estimates  was  used  to  yield  the  final  decision.  In 
this  chapter,  we  still  consider  the  problem  of  multisensor 
CFAR  detection  with  data  fusion.  However,  the  compressed 
data  transmitted  between  sensors  are  decisions  instead  of 
estimates.  In  Section  6.2,  we  formulate  the  general  problem 
of  adaptive  CFAR  detection  for  two  different  network  topolo¬ 
gies.  These  systems  are  optimized  and  their  performance  is 
studied.  In  Section  6.3,  we  illustrate  the  performance  by 
an  example.  In  Section  6.4,  we  provide  a  summary  along  with 
a  discussion. 


The  network  topologies  to  be  considered  in  this  chapter 


are  shown  in  Figures  6-1  and  6-2.  In  the  first,  decisions  from 
individual  CA-CFAR  processors  are  transmitted  in  parallel  to 
the  CFAR  processor  where  a  final  decision  is  to  be  made.  In 
the  second,  the  CFAR  processors  are  connected  in  tandem.  De¬ 
cisions  are  transmitted  from  one  CFAR  processor  to  the  next. 

At  the  last  CFAR  processor,  the  final  decision  is  made. 

The  target  to  be  detected  is  assumed  to  be  of  Swerling 
type  I  embedded  in  a  white  Gaussian  noise  of  unknown  level. 

At  each  CA-CFAR  detector,  CDi,  i  =  l,2,...,n,  we  denote 
the  target  SNR  by  S • ,  i=l,2,...,n.  The  target  SNR  at  the 
final  CA-CFAR  detector,  CDq,  is  denoted  by  S.  For  simplicity 
and  without  loss  of  generality,  we  let  S-^  =  S2= .  .  .  =Sn=S .  The 
number  of  range  cells  used  by  each  CA-CFAR  detector,  CD^  is 
denoted  by  N^,  i=l,2,...,n,  and  N  denotes  the  number  of 

range  cells  for  detector  CD  .  Next,  we  consider  the  two 

o 

networks  topologies. 


6.2-1.  Adaptive  CFAR  Detection  with  Parallel 
Network  Topology 


The  system  under  consideration  is  shown  in  Figure  6-1 
It  consists  of  (n+1)  CA-CFAR  processors.  Based  on  its  ob¬ 


servation,  each  CA-CFAR  detector,  CD^,  i=l,2,...,n,  makes  a 


decision  D^,  i=l,2,...,n,  where  D^  c  {0,1}.  The  adaptive 
threshold  of  CD^  is  scaled  by  the  threshold  multiplier 
T.,  i=l,2,...,n.  As  will  be  shown  later,  these  T,  rs  are 


Tandem  Network  Topology  tor  Adaptive  ltak 
Detection . 


obtained  so  as  to  maximize  the  overall  probability  of  detec¬ 
tion  while  maintaining  the  overall  probability  of  false-alarm 
at  the  desired  value.  The  CA-CFAR  detector,  CD  ,  combines 
the  decisions  D^,  i=l,2,...,n,  with  its  own  estimate  of  the 
level  of  the  background  noise  to  obtain  the  adaptive  threshold 
Q.  The  threshold  Q  is  scaled  by  a  threshold  multiplier,  which 
is  also  obtained  so  as  to  maximize  the  overall  probability  of 
detection  while  CFAR  is  achieved.  The  threshold  multiplier 
depends  on  the  decisions  D^,  i=l,2,...,n.  We  define  the 


vector,  DV^,  consisting  of  decisions  arriving  at  CDq  as 


DV^  =  (D1,D2, . . . ,Dn)T  ,  £  =  l,2,...,2n 


(6-1) 


Since  each  ,  i  =  l,2,...,n,  takes  the  value  zero  or  one, 

DV^  can  take  2n  possible  values.  For  each  value  of  DVjj, 
arriving  at  CDo ,  there  is  a  corresponding  threshold  multiplier 
TqJI,  £  =  l,2,...,2n,  at  CDq  .  These  threshold  multipliers, 

Tq£,  £  =  l,2,...,2n,  are  obtained  by  optimizing  the  overall 
system,  i.e.,  by  maximizing  the  overall  probability  of  detec¬ 
tion,  while  the  overall  probability  of  false-alarm  is  kept 
constant.  Thus,  the  threshold  at  CDq  adapts  to  two  types  of 
information,  its  own  estimate  of  the  level  of  the  background 

noise  and  the  incoming  decision  vector  DV„ .  Let  Pn  ,  P_  ,  and 

^  i  i 

P^  be  the  probability  of  detection,  the  probability  of  false- 
1 

alarm, and  the  probability  of  a  miss  at  each  CA-CFAR  detector, 
CD^,  i=l,2,...,n,  respectively.  The  overa 1 1  probability  of  detec¬ 
tion,  the  overall  probability  of  a  miss,  and  the  overall  probability 
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of  false-alarm  at  the  output  of  the  CA-CFAR  detector,  CDq, 
are  denoted  by  P^,  P^,  and  Pp,  respectively.  Then,  PD  and  Pp 
can  be  written  as 


Pn  -  l  P„  n_  PM  n,  Pn 
D  i-1  D  S°  Mj  S1  Dk 


Pp  ■  l  PF  "o  (1  •  PF  >  "l  PF 

h  t=l  b  S°  j  s1  hk 


(6-2) 


(6-3) 


where  from  Appendix  A, 


pd  ’  PrCDo  ‘  Hi> 


r  1+S  ,N 


PF  '  Pr(Do  '  UOV  Hc> 


(1  +  To*> 


(6- 4 . a) 


(6-4. b) 


(6-5. a) 


(6-5. b) 


•  £  9 

The  superscript,  l,  in  p^  and  Pp  indicates  the  correspondence 
between  the  probabilities  of  detection  and  faise-alarm  and 
the  specific  threshold  multiplier,  Tq^,  obtained  from  the 
optimized  system. 

To  optimize  the  overall  system,  we  form  the  following 
objective  function. 


J  “  PD  +  ^PF 


(6-6) 
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Substituting  equations  (6-2),  (6-3),  (6-4. b)  and  (6-5. b) 


into  equation  (6-6),  we  get 
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^  1 
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1=1  (i+t olr 
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N. 

(1+TM?)  1 


N. 

(1  +  S)  1 


1  *  nl  - ~ -^-N. 

s  (1+S+TmJ)  1 


(1+TMt  )  1 


(6-7) 


where  TM.'s  are  the  threshold  multipliers  of  CD.'s, 
i=l,2,...,n.  The  superscript  zero  in  TM°  indicates  that  the 
decision  of  the  corresponding  CD^  is  zero,  and  the  super¬ 
script  one  in  TM^  indicates  that  the  decision  of  the  cor¬ 
responding  CD^  is  a  one.  In  order  to  maximize  Pp  while  Pp 
is  fixed,  for  each  value  of  DV^,  we  must  find  the  best 
threshold  multiplier,  T  at  CDq  and  the  best  set  of  thresh 
old  multipliers  at  the  CD^ ' s ,  i=l,2,...,n.  Maximizing  Pp 
for  a  fixed  Pp  is  equivalent  to  maximizing  the  objective 
function  J.  Thus,  we  take  the  derivatives  of  J  with  respect 
to  TM°  and  TM^,  9.  =  1 , 2  ,  .  .  .  ,  2n ,  i=l,2,...,n  (the  number  of 
TM° ' s  and  TM^'s  add  to  n) ,  and  set  them  equal  to  zero.  The 
resulting  equations  are 
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6.2-2.  Adaptive  CFAR  Detection  with  Tandem 
Network  Topology 

In  this  case,  the  system  considered  is  shown  in 
Figure  6-2.  The  CA-CFAR  detector,  CD^ ,  makes  a  decision, 

D^,  and  transmits  it  to  the  CA-CFAR  detector  CE>2'  Based 
on  the  decision  and  on  its  own  estimate  of  the  background 
noise,  CA-CFAR  detector  CD^  makes  another  decision,  D^.  The 
decision  D2  is  transmitted  to  the  next  CA-CFAR  detector  CD^, 
and  the  process  continues.  At  the  last  CA-CFAR  detector, 

CDq,  the  final  decision,  Dq ,  is  made. 

Let  the  threshold  multipliers  of  detectors  CDq,  CD^, 

CD2»...,CD  ,  be  denoted  by  T,  T^,  T2»...,T  ,  respectively. 

We  define  P  ,  P  ,  and  P  ,  i=l,2,...,n,  to  be  the  proba- 
ui  hi  Mi 

bility  of  detection,  the  probability  of  false-alarm,  and 
the  probability  of  a  miss  at  the  individual  CA-CFAR  detectors 
CD^,  i=l,2,...,n.  At  CD^,  i=2,3,...,n,  the  threshold  multiplier 
is  represented  by  TV(Tl)  if  the  decision  received  from  CD.^ 
is  zero  (one).  Let  N,  N^,N2,...,Nn  be  the  number  of  range 
cells  used  by  CA-CFAR  detectors  CDq ,  CD^ ,  CD2 , . . . , CD^ ,  respec¬ 
tively.  As  before,  we  assume  that  the  target  to  be  detected 
is  of  Swerling  type  I  and  has  an  SNR  S  for  all  of  the  detec¬ 
tors.  Then,  the  probability  of  detection  and  the  probability 


of  false-alarm  of  CA-CFAR  detector  CD^  are 


(1  +  S) 


(6-11) 


(1+S+TX) 


1  1 
1  (1+TX)  1 


(6-12) 
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The  probability  of  detection  and  the  probability  of  false- 
alarm  at  CA-CFAR  detector  CD2  are 


PD  =  Pr(D2=l |D1=1,  H:)  PrCD1=l|H1)  + 
Pr(D2=l|D1=0,  Hx)  Pr(D1»0|H1) 


(6-13. a) 


Pr(D2  =  l|D1  =  l>  H^P^  +  Pr(D2  =  l|D1=0,  (1-P^) 


(6-13. b) 


Pp  =  Pr(D2=l|D1=l,  Hq)  Pr(D1=l|Ho)  + 
Pr(D2=l|D1=0,  Hq)  Pr(D1=0|Ho) 


(6-14. a) 


=  Pr(D2  =  l|D1=l,  Ho)Ppi  ♦  Pr(D2=l|D1  =  0,  Hq)  (1-P^) 


(6-14. b) 


The  resulting  expressions  for  P^  and  Pp^  are 


>  =  - - py—  +  —  u raj — _  n - ■ 

D  j  N?  NT  NT  fF 

(1+S+T2)  z(1+S+T1)  1  (1+S+T°)  L  (1+S+Tj )  1 


Nl+N2 


*2  (1+T2)”2(1+T1)”1  (1+T°)  2 

For  CA-CFAR  detector  CD3>  we  have 


(6-15) 


-]  (6-16) 


(1+T1) 


PD3  =  Pr(D3=l|D2=l,  Hx)  Pr(D2=l)H1)  + 

Pr(D3=l|D2=0,  Hj)  Pr(D2=0|H1)  (6-17. a) 
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PF3  =  Pr(D3=l|D2=l,  Ho)  Pr(D2=l|Ho)  + 

Pr(D3=l |D2=0,  Hq)  Pr(D2=0|HQ)  (6-18. a) 


— Vfq  rK  *  —Vic  (1-pf2) 

(1+T3)  Z  (1+T°)  1 


(6-18. b) 


where,  P*  =  Pr(D2  =  1  |  ,  P°  =  Pr  (D2  =  0  |  Dj_ ,  , 

P1  •  Pr(52  -  1  IDj,  Ho)  and  P°  -  Pr(B2  =  0^,  Hq).  pJ 

and  P  are  equal  to  equations  (6-15)  and  (6-16)  respectively, 
2  o  o 

while  P°  and  P°  can  be  expressed  as 
2  2 


. a  +  s)  ,  p  + 

PV  D1 

(1+S+T2)  z  1 
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(1+S+T°)  c  1 

(6-19) 


p“f2  ■  I1  -  pFl  *  t1  -  —7^]  n-pF  1 

(1  +  K)  1  (1+T°)  2  1 

Z  (6-20) 

The  superscript  one  (zero)  in  pi  (P°  )  and  P^  (P°  ) 

U2  U2  F2  F2 

indicates  that  the  decision  made  by  detector  two  is  one  (zero). 
Recall  that  the  superscript  one  (zero)  in  T2  (T°)  indicates 
that  the  received  decision  from  is  one  (zero).  Continuing 
in  this  manner,  we  can  write  the  expressions  of  the  probability 
of  detection  and  the  probability  of  false-alarm  at  the  ith 
detector  as, 


IT1  pFl  +  t1  ' 


=  Pr(Di  =  l|D._1  =  l,  Hp  Pr(Di_1=l  |H1)  + 
Pr(Di=l|D._1=0,  Hx)  Pr(Di_1-0|H1) 


(6-21. a) 
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Pp  =  Pr(Di=l |Di_1*l,  Hq)  Pr(D._1=l|Ho)  ♦ 


Pr(Di=l|D._1=0>  Hq)  Pr(Di.1=0|Ho) 


1  pi  t  1  ( ■)  .  pQ  ) 

- nr  pf .  .  ~r  f-  i 
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where 
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,  j-0,1. 


(6- 22. c) 
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,  j  =  0, 1 


(6-22. d) 


The  superscript  j  in  P^.  and  P^  indicates  whether  the 

i-l  hi-l 

previous  decision  is  zero  or  one.  The  general  expressions 
for  the  overall  probability  of  detection,  Pp,  and  the  overall 
probability  of  false-alarm,  Pp,  at  the  CA-CFAR  detector  CDq 


can  then  be  written  as 


P  =  (1+S)N  pl  +  (1+S)N 

D  (l+S+T1)^  Dn  (1+S+T°)f 


(1-P?  > 
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(6-23) 


PF 


(l+T1) 


— Vn  ci-Pp  ) 

(1+T°)N  Fn 


(6-24) 


where,  Pp  and  P^  ,  j=0,l,  are  defined  recursively  in  equa- 
n  n  in 

tion  (6r22.c)  and  (6-22. d)  and  T  (T  )  is  the  threshold  at 

detector  CD  if  the  received  decision  D  is  one  (zero), 
o  n 


We  wish  to  maximize  the  probability  of  detection,  P^,  at 
the  output  of  CA-CFAR  detector  CDq,  while  the  probability  of 
false-alarm,  Pp,  is  fixed  at  the  desired  value  v.  Since  the 
CA-CFAR  detectors  CDq,  CD^,  CD2»...,CD  are  coupled,  the  en¬ 
tire  system  is  optimized.  To  achieve  this,  we  form  the 
objective  function 


J  -  PD  *  5  [pF 


(6-25) 


We  take  the  derivatives  of  J  with  respect  to  T-*  and  T^, 
j=0,l,  and  k=l,2,...,n,  and  set  them  equal  to  zero.  This 
results  in  a  system  of  2(n+l)  nonlinear  equations  in  2fn+l) 
unknowns  to  be  solved.  In  the  next  section,  we  present  an 
example  to  illustrate  the  effectiveness  of  the  system  under 
consideration. 


6.3.  Example 


The  system  under  consideration  is  shown  in  Figure  6-3. 
In  order  to  avoid  any  ambiguity  with  the  notation  defined 
in  the  subsections  6.2-1  and  6.2-2,  we  define  a  more  explicit 
notation  that  is  used  in  this  example. 


D  ,D-  =  detection  decisions  at  detectors  CD  and 
o  1  o 

CDp,  respectively. 

Tp  =  threshold  multiplier  at  CA-CFAR  detector  CD^ 
T°  =  threshold  multiplier  at  CA-CFAR  detector  CD 


Phenomenon  H 


Fig.  6-3.  A  Two  Sensor  Network  Topology 
for  Adaptive  CFAR  Detection. 
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T1  m  threshold  multiplier  at  CA-CFAR  detector 

CD  if  the  decision  D.  is  one. 
o  1 

^D’^D  =  Pr°t,a^:*-^ities  °f  detection  at  CA-CFAR 

detectors  CDq  and  CD^,  respectively. 

Pp.Pp  =  probabilities  of  false-alarm  at  CA-CFAR 
F  F1 

detectors  CDq  and  CD^,  respectively. 

S  =  target  SNR  at  CA-CFAR  detectors  CDq 
and  CD^ 

N,N^  =  number  of  reference  cells  at  CA-CFAR 
detectors  CDq  and  CD^,  respectively. 

PD  and  Pp  can  then  be  expressed  as 


PD  =  Pr(Do=l|D1=l,  Hx)  PrCDp-llHj)  + 
Pr(Do=l|D1=0,  Hx)  PrCDj-OlHp) 

PF  =  Pr CDq=1 |D1=1,  Hq)  Pr(D1=l|Ho)  + 

PrCDo=l|D1=0,  Ho)  Pr(D1=0|Ho) 


(6-26) 


(6-27) 


Substituting  the  expressions  of  Pr (Dq= 1 | D^  =  l ,  H^)  , 

Pr(Do=l|D1  =  l,  Ho),  Pr(Do=l|D1  =  0,  ,  Pr  (DQ«  1 1  D^O ,  Hq)  , 

Pn  ,  P  ,  P  and  Pn  into  equations  (6-26)  and  (6-27), 
L»1  r-^  ^2  u2 

P„  and  P„  become 
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(1  +  TX) 
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Substituting  equations  (6-28)  and  (6-29),  and  Pp  -  v  into 
equation  (6-25),  the  objective  function  is 


(1+S+T1)N(1+S+T1)  1  (1+S+T  ) 


( 1+S+Tp ) 


+  ? { — -~r - TT  +  [1  ~  nt 

(l+T1)N(l+T1)  1  ^  T  ->  (1  +  Tp)  1 


]-v}  (6-30 


Maximizing  Pn  is  equivalent  to  maximizing  J.  That  is,  we 


-1  0  rri  1 


take  the  derivatives  of  J(T  ,  T  ,  T^j  with  respect  to  T  , 
and  T, ,  respectively,  and  set  them  equal  to  zero  to  obtain 


(1+S+T  ) 


(1+S+T, ) 


(1+S+T1)N+1 (1+S+T, )  1  (1+TX)N+ 


The  equations  (6-31),  (6-32),  (6-33)  and  Pp  =  v  form  a  set 
of  four  nonlinear  equations  in  four  unknowns  to  be  solved. 
We  first  obtain  the  Lagrange  multiplier,  from  equation 


ss 

of  detection,  ,  of  the  single  sensor  system,  for  the  same 
probability  of  false-alarm  and  the  same  target  SNR. 


CA-CFAR  detectors  whose  operation  is  as  described  in  this 
chapter.  For  the  same  cross  sectional  target  SNR  S,  for  an 
overall  probability  of  false-alarm  Pp  =  10  \  for  N^  =  6  and 
N2=4,  we  compare  the  performance  of  the  three  networks  in 
Table  6~2. 


s 

Pa 

Pb 

D 

PC 

D 

5 

0.0152 

0.0152 

0.0687 

10 

0.166 

0.166 

0.179 

15 

0.  526 

0.  526 

0.594 

Table  6-2.  CA-CFAR  Performance  of 
the  Two-Sensor  Systems 
of  Figure  6.4. 

The  superscripts  a,  b  and  c  on  the  probabilities  of  detection 
of  the  entries  in  Table  6-2,  represent  the  systems  a,  b, 
and  c  as  shown  in  Figure  6-4,  respectively.  It  should  be 
noted  that  if  we  had  used  the  "OR"  fusion  rule  instead  of 
the  "AND"  fusion  rule  in  the  system  shown  in  Figure  6-4(a), 
the  performance  in  terms  of  P^  would  have  been  worse  for  the 
values  of  S  considered  here. 

6.5.  Summary  and  Conclusions 

In  this  chapter,  we  have  derived  expressions  for  the 
probability  of  detection  using  n  CA-CFAR  detectors  in  parallel 


1  58 


‘t, 


and  in  tandem.  In  the  parallel  topology,  n  decisions  from 
n  CA-CFAR  detectors  arrived  at  the  (n+l)th  CA-CFAR  where  the 
final  decision  was  obtained.  The  entire  system  was  optimized 
and  the  best  threshold  multipliers  at  the  CA-CFAR  detectors 
were  obtained  so  that  the  overall  probability  of  detection 


was  maximum  while  the  overall  probability  of  false-alarm  was 
maintained  at  the  desired  value.  In  the  tandem  topology,  the 
detection  decisions  were  transmitted  from  one  CA-CFAR  detector 
to  the  next  in  series  and  the  entire  system  was  optimized. 

The  improvement  in  the  performance  was  shown  by  an  example. 
This  procedure,  in  which  the  threshold  is  adaptive  to  two 
types  of  compressed  data,  the  received  decision  and  its  own 
estimate  of  the  level  of  the  background  noise,  turned  out  to 
be  very  effective. 


139 


hji  nji  r-jL  TVjf  r.M  w  JC  AJi  AX FJ*  7 


r^jnc^^r 


CHAPTER  VII 


SUMMARY  AND  SUGGESTIONS  FOR  FUTURE 


RESEARCH  WORK 


7.1.  Summary 

New  adaptive  threshold  techniques  for  CFAR  processing 
were  proposed  and  analyzed.  In  the  first  half  of  the 
report,  we  studied  CFAR  processors  for  multiple  target  situ¬ 
ations.  We  proposed  a  weighted  cell-averaging  CFAR  detector 
for  multiple  target  situations.  We  obtained  the  optimum 
weights  for  the  leading  and  the  lagging  reference  windows  so 
that  the  probability  of  detection  was  maximum  while  the 
probability  of  false-alarm  was  kept  at  the  desired  value. 

The  performance  of  the  WCA-CFAR  detector  was  shown  to  be 
superior  to  that  of  the  CA-CFAR,  GO-CFAR  and  SO-CFAR  de¬ 
tectors.  The  WCA-CFAR  detector  was  shown  to  be  especially 
effective  when  the  interference  is  large. 

The  second  CFAR  detector  proposed  for  a  multiple 
target  environment  was  a  ce 1 1  -  censored  CFAR  detector.  A 
predetermined  fixed  threshold  at  the  level  of  each  range 
cell  was  used  to  determine  the  presence  or  absence  of  inter¬ 
ference  in  each  range  cell.  The  cells,  for  which  inter¬ 
ference  was  decided,  were  censored  while  others  (cells 
without  interference)  were  used  in  the  estimation  of  the 
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background  noise.  The  sum  of  the  leading  cells  without 
interference,  and  the  sum  of  the  lagging  cells  without  inter¬ 
ference  were  added  to  yield  the  adaptive  threshold,  which  in 
turn  was  scaled  to  achieve  CFAR.  The  performance  was  studied 
and  it  was  shown  that  the  cell-censored  CFAR  detector  was 
effective  in  a  multiple  target  environment. 

In  the  second  part  of  the  report,  we  developed 
the  theory  of  distributed  adaptive  CFAR  detection.  CFAR  de¬ 
tection  with  distributed  detectors  and  data  fusion  was 
studied.  Given  a  "k  out  of  n"  fusion  rule  at  the  data  fusion 
center,  the  optimum  threshold  multipliers  at  the  CFAR  detec¬ 
tors  were  obtained  so  that  the  overall  probability  of  detec¬ 
tion  at  the  output  of  the  data  fusion  center  was  maximum, 
while  the  overall  probability  of  false-alarm  was  fixed.  Then, 
the  entire  system  was  optimized,  i.e.,  the  optimum  threshold 
multipliers  at  the  individual  detectors  were  obtained  as  well 
as  the  optimum  fusion  rule  at  the  data  fusion  center.  In  the 
optimization  of  the  entire  system,  the  overall  probability  of 
false-alarm  was  maintained  at  the  desired  value  while  the 
overall  probability  of  detection  was  maximized  at  the  data 
fusion  center. 

Due  to  the  fact  that  more  knowledge  about  the  back¬ 
ground  noise  improves  the  detection,  adaptive  CFAR  detection 
with  multiple  background  estimators  was  studied.  Estimates 
of  the  background  noise  were  transmitted  from  the  estimators 
to  the  CFAR  detector.  The  CFAR  detector  computed  its  own 
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estimate  of  the  background  noise  and  combined  it  with  the 
received  estimates  to  yield  the  adaptive  threshold.  The 
adaptive  threshold  was  then  scaled  and  was  compared  with  the 
statistic  of  the  output  of  the  cell  under  test  to  make  a  de¬ 
cision.  As  expected,  the  performance  improved.  An  exten¬ 
sion  using  a  hybrid  system  consisting  of  background  esti¬ 
mators  and  a  data  fusion  center  was  also  considered.  The 
performance  of  the  hybrid  system  was  also  studied  by  means 
of  an  example. 

Finally,  adaptive  CFAR  detection  for  two  different 
network  topologies  was  considered.  The  two  topologies  con¬ 
sidered  were  the  parallel  and  the  tandem  topologies.  In 
these  cases,  the  compressed  data  that  were  transmitted  from 
one  CFAR  detector  to  the  next  were  decisions  instead  of  esti¬ 
mates.  The  final  decision  at  the  last  CFAR  detector  was 
based  on  two  types  of  compressed  data;  namely,  its  own  esti¬ 
mate  of  the  background  noise  and  the  received  decisions.  In 
both  cases,  the  entire  system  was  optimized  so  that  the 
probability  of  detection  was  maximum,  while  the  probability 
of  false-alarm  was  maintained  at  the  desired  values.  In 
maximizing  the  probability  of  detection,  optimum  threshold 
multipliers  were  derived.  Again,  as  expected  the  perfor¬ 
mance  improved.  In  the  next  section,  we  provide  some  sug¬ 
gestions  for  future  research  work. 


Suggestions  for  Future  Research  Work 


As  indicated  earlier,  no  work  has  been  reported  in  the 
literature  on  CFAR  detection  using  distributed  sensors.  In 
this  report,  we  make  the  first  contribution  in  this  area 
by  developing  the  theory  of  distributed  adaptive  CFAR  detec¬ 
tion.  We  did  not  consider  multiple  target  situations  in  our 
work.  One  obvious  direction  for  future  research  is  to  extend 
the  theory  developed  here  for  multiple  target  situations. 
Also,  different  network  topologies,  structures  and  implementa 
tions  can  be  considered  and  analyzed  for  distributed  adaptive 
CFAR  detection. 

Also,  as  shown  in  Chapters  four  and  six,  solving  the 
nonlinear  equations,  to  obtain  the  optimum  threshold  multi¬ 
pliers  at  the  individual  CFAR  detectors  and/or  the  optimum 
fusion  rule  at  the  data  fusion  center,  required  extensive 
computer  iterative  steps.  These  values  may  be  computed  and 
compiled  in  tables  so  that  a  table  look-up  procedure  may  be 
implemented.  Development  of  more  efficient  numerical  tech¬ 
niques  for  solving  the  above  nonlinear  equations  is  desir- 
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APPENDIX  A 


A- 1 .  Swerling  Cross  Section  Models  [12,36] 

Case  1.  In  this  case,  the  returned  signal  power  per  pulse 
on  any  one  scan  is  assumed  to  be  constant,  but  these  echo 
pulses  are  independent  (uncorrelated)  from  scan  to  scan. 
This  assumption  ignores  factors  such  as  the  effect  of  the 
antenna  beam  shape  on  the  amplitude  of  the  returned  signal. 
A  returned  signal  of  this  type  is  referred  to  as  scan-to- 
scan  fluctuation.  The  probability  density  function  for  the 
cross  section  (input  signal -to-noise  power  ratio,  a,  is 


assumed  to  be 


pCa>aav)  =  T~  exp  C"  ’  a  -  0 

av  av 


(A- 1 ) 


where  oav  is  the  average  cross  section  (average  of  a)  over 
all  target  fluctuations. 


Case  2.  In  this  case,  the  fluctuations  are  more  rapid  than 
in  case  1  and  are  assumed  to  be  independent  from  pulse  to 
pulse  instead  of  from  scan  to  scan.  This  type  of  fluctu¬ 
ation  is  referred  to  as  pulse-to-pulse  fluctuation.  The 
probability  density  function  for  the  target  cross  section  is 
the  same  as  given  in  equation  (A-l). 
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Case  3.  In  this  case,  the  fluctuation  is  as  in  Case  1, 
i.e.  scan-to-scan  fluctuation,  but  the  probability  density 


function  is  given  by 


P  (a  ,  a  )  = 
v  '  av' 


2a/ a. 


,  a  >  0 


(A-  2) 


Case  4.  In  this  case,  the  fluctuation  is  pulse-to-pulse  as 
in  case  2,  but  the  probability  density  function  is  given  by 
equation  (A- 2). 

Only  Case  1  is  considered  in  this  dissertation. 


A- 2 .  Reference  Window  Statistic  [37] 

When  the  background  noise  is  white  Gaussian  of  un¬ 
known  level,  the  probability  density  function  of  the  out¬ 
put  at  each  range  cell  is 


1 

Pq  Cq±)  =  j  e  »  ;!  ^ 


(A-3) 


- 7  ^  w  +>  j 

^i  1  2a^  1 

2 

where  a  is  the  average  noise  power.  In  cell-averaging,  the 
outputs  of  these  range  cells  are  assumed  to  be  statistically 
independent.  A  reference  window  U  of  N  range  cells  is  the 
sum  of  all  these  cells,  i.e., 


U  =  l  q. 
i  =  l  1 


(A-4) 


where  Pn  (q.),  i=l ,2 , . . . ,N, is  as  defined  in  equation  (A-3) 

Vi  1 

To  find  Pjj(n)  we  use  the  following  theorem  from  [37]. 
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rheorem.  If  q^ ,q2 , • • • ,q^  are  mutually  independent  random 
variables  with  the  exponential  distribution  function  as  in 
equation  ( A -  3 ) ,  then  the  sum  U  =  +...+  q^  has  a  proba¬ 

bility  density  function  Pjjfu)  and  a  distribution  function 
Fu(u)  given  by 


Vu) 


FyCu)  =  1-e 


1  r  u  -,N-1  -u/2 a  ^  n 

rroy  e  ’  u  - 0 


CA-5) 


,,  2  ,  ,  N-l 

u/2°  (1  +  u_  t  (  u  ,2  .  1  u  , 

2o  *1!  2a  c'  2o 


1 _ ] 

(N-l)  ) 


u  >  0 


(A- 6) 


Proof .  We  prove  this  theorem  by  induction.  For  N=1 ,  the 
theorem  holds  because  equation  (A-5)  simply  reduces  to  equa¬ 
tion  (A-4).  Now,  assuming  that  the  theorem  is  true  for  N, 
we  will  show  that  it  is  also  true  for  N+l.  The  density 
Py(u)  for  U  defined  in  equation  (A-4)  is  given  by  the  follow¬ 
ing  convolution 


wmt 

ME. 


Pfl+lCt)  =  j  PNCt-q)  Pj(u)  du 


(A-  7) 


which  reduces  to 


_  (l/2a2)N+1  -t/2 a2 
(N-l) !  e 


1  (t/2o2)N  - t/2a2 

TT  N! 

2a 


CA-8) 


quation  (A-8)  is  of  the  form  of  equation  (A-5)  which  means 
tat  equation  (A-5)  holds  by  induction  for  all  N.  Taking 
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the  derivative  of  equation  (A-6) ,  we  obtain  equation  (A-5) 
which  justifies  that  equation  (A- 6)  holds. 


A-3 .  Detection  Probability  Without  Interference  in 
the  Range  Cells 

In  the  conventional  CA-CFAR  with  N  reference  cells  and 
reference  windows  U,  V,  the  normalized  probability  density 
functions  of  U  and  V  are 


-  1 

V»)  =  — TT  “  e‘U 

u  r(|) 


e  ,  u  >  0 


*  i 

n  ,  ,  1  7 "  1  -v 

Vv)  'rx1 

1  c7} 


e  ,  v  _>  0 


(A-9) 
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The  threshold  Q  is  the  sum  of  U  and  V,  and  therefore,  the 
probability  density  function  of  Q  is  given  by  the  convolu¬ 


tion 


Pq(q)  =  PyCvij  Pv(q-u)  du 


n  q  N  .  n  - 

e-q  r  7-1 

- n — 2  I  u  Cq-u)^  du 

fr  C7)  1 2  i 


lr (^) 1 2 


e  q  qN" 1  bCy*  7) 
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The  probability  of  detection  for  a  target  SNR  S  of  Swerling 


Case  1  is 


"  -  Tq 

pd  =  e  1+s  PqCq)  dq 
0 


(A- 1 2) 


From  [38],  we  obtain 


p  =  r _ _ _iN 

*D  ll  +  S  +  TV 


(A-13) 


A- 4.  A  Reference  Window  Statistic  With  One 
Interfering  Target  [14] 

If  the  reference  window  is  U  with  one  interfering 
target,  having  an  SNR  I,  of  Swerling  type  I  and  N  range 
cells,  U  can  be  written  as 


where , 


N-l 

u  -  l  *  qN 


p0  (qj)  ■  e  1,  qj  >  o  , 


i=l , 2 , . . . ,N-1 


(A-14) 


PQ*,CqN}  =  1  +  1  6 


— 


(A-14 . a) 


(A-14 .b) 


Note  that  we  have  assumed  that  the  interference  is  in  the 
Nth  cell,  but  it  could  be  in  anyone  of  the  N  range  cells. 
Let , 
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then,  from  equation  (A-2) 


PX(X)  =  (N-2TT  xN"2  e"X  ’  X-°  (A‘17^ 


Since  U  is  the  sum  of  X  and  Y,  the  probability  density  func¬ 


tion  of  U  is 


Py(u)  =  j  pxCx)  py (u"x)  dx 

0 


(A-18) 


Substituting  equations  (A-14.b),  (A-16)  and  (A-17)  into 
equation  (A-18) ,  we  obtain 


PuCu)  = 


1  N-2  -x  - (u-x)/l+I  . 

•  x  e  e  dx 
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From  [38]  ,  we  have 


n  -  viXj 
x  e  dx 
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[u  >  0,  Re  v  >  0] 


(A- 19) 


(A- 20) 
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Thus,  using  equation  (A-20)  inequation  (A-19)  and  re¬ 
arranging  terms,  the  probability  density  function  of  U  is 


pu(u) 


1  J-1  +  I-.N-2 
I  e 


•u  Nr2  1  r  I  ,j  j 

Jo  IT  (TTTlJ  u  • 


U  >  0 


(A- 21) 


A-5.  Some  Values  of  a  and  6 


The  values  of  a  and  S  are  obtained  by  solving  the 
nonlinear  equations  (2-43),  (2-44)  and  (2-46).  The  values  of 
a  and  3  for  the  cases  considered  are  given  below. 


Pr  =  10 
r 


I  =  S  x  10,  N=4 


s 

a 

S 

5 

1.4318 

8.4275 

10 

.2711 

15.738 

15 

.0680 

19.412 

20 

.0207 

20.233 

25 

.  0065 

20.451 

30 

.0021 

20.515 
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MISSION 

of 

Rome  Air  Development  Center 

KAVC  plan*  and  executes  research,  development,  test 
and  selected  acquisition  p4.ogA.ams  In  suppoat  of 
Command,  Control,  Communications  and  Intelligence 
(C^I)  activities.  Technical  and  engineering 
support  within  areas  o f  competence  Is  provided  to 
ESV  Program  Offices  (POs)  and  other  ESV  elements 
to  perform  elective  acquisition  of  C^I  systems. 

The  areas  of  technical  competence  Include 
communications ,  command  and  control,  battle 
management,  Information  processing,  surveillance 
‘ sensors ,  Intelligence  data  collection  and  handling, 
solid  state  sciences,  electromagnetics ,  and 
propagation,  and  electronic,  maintainability , 
and  compatibility. 
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